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Abstract 
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■ The charge and magnetic form factors of hght nuclei, mainly for mass number 

, A<4, provide a sensitive test of our understanding of nuclei. A number of "exact" 

calculations of the wave functions starting from the nucleon-nucleon interaction are 
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' available. The treatment of two-body effects needed in the calculation of the elec- 

. tromagnetic form factors has made significant progress. Many electron scattering 

I experiments have provided an extensive data base from which the various (mainly 

^ ■ elastic) form factors can be extracted. This review discusses the data and the de- 
termination of the form factors, and compares them to the results of theory. 
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^ '. 1 Introduction 

Few-body nuclei have long attracted particular attention in nuclear physics. This essen- 
tially is due to the fact that for these nuclei non-relativistic and sometimes relativistic 
calculations of the wave function with high accuracy starting from the nucleon-nucleon 
(N-N) interaction can be performed. For heavier nuclei models, or approaches with various 
approximations, have to be employed as a consequence of the complexity of many-body 
systems. In this case one largely looses the connection to the underlying, microscopic, 
physics. 

For these light nuclei one then can best test our understanding of nuclei. Can they be 
described as essentially non-relativistic systems of nucleons interacting via the known N-N 
force, or do we have to go beyond this picture which has become known as the "standard 
model of nuclear physics" ? To which degree are nuclei relativistic systems? Do we need to 
explicitly allow for other constituents like excited nucleons, pions,..? Does the composite 
(quark) structure of the nucleon play an important role? 

The few-body nuclei provide the possibility to study this rather wide spectrum of 
questions. While the deuteron is a very dilute system, the '^He nucleus is the nucleus 
with largest density. In contrast to heavy nuclei, these nuclei have large gradients of the 
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density which also can play an interesting role. The A=3 nuclei offer the possibility to 
study near-identical systems that differ only by the exchange of the role of protons and 
neutrons. 

Significant progress has been made during the last decade in quantitatively computing 
the wave function of light nuclei. For many observables one today can be sure that 
differences in the calculated results reflect genuine differences in the physics included, and 
not simply approximations that needed to be made in the course of the calculation. A 
good part of these advances obviously is due to the vastly increased computational means 
that have become available. 

Electron scattering provides an excellent tool for a detailed check of theoretical calcu- 
lations of the wave functions of few-body nuclei, a test which often is more informative 
than what can be learned from integral observables such as energies of states or moments. 
Provided the experimental (e, e) data go to large (four) momentum transfer q, the spatial 
structure of the one-body densities can (in Impulse Approximation lA) be studied with 
good spatial resolution, of the order of 1.5/g. 

As a consequence of the great interest in few-body form factors, many experiments 
on few-body nuclei have been performed over the last decades despite the experimental 
challenge they represent. We today have a rather complete set of data to which one can 
compare. An important aspect of the recent progress in (e, e) experiments is the ability 
to measure polarization observables in the region of q where they can help to discriminate 
between different theoretical predictions. Much of this progress has been made possible 
by the modern high-energy electron accelerators with high duty cycle such as MAMI or 
JLAB and the development of polarized sources, targets and polarimeters. 

The form factors of light nuclei are of particular interest for another reason: Selected 
form factors, e.g. the isovector magnetic ones, exhibit a special sensitivity to exchange 
current contributions. For heavier nuclei the (magnetic) form factors show diffraction 
features at much lower momentum transfer already due to the structure of the wave 
function, and the contribution of two-body currents is much more difficult to isolate. 

There is quite a vast literature on the few-body form factors — too vast for an ex- 
haustive list of references without omitting significant work. These reviews have different 
areas of emphasis. The present review — written by an experimentalist — gives more 
than most others attention to data and ignores much of the detailed formalism of prime 
interest to more theory-oriented reviews. For general reviews on electron scattering see 

1 1,1- 

In this review, we first discuss the ingredients necessary for the calculation of the form 
factors. We then discuss for A=2,3,4,6,16 the data available, the determination of the 
optimal set of experimental form factors and their comparison to theoretical calculations. 
For the A=2 case, we include the transition form factor to the ^S state. 

2 Ingredients 
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2.1 Nucleon form factors 



The charge- and magnetic form factors of nuclei are calculated from the wave function of 
the nucleons, obtained in most cases from a solution of the Schrodinger or Bethe-Salpeter 
equation for pointlike nucleons. In order to obtain the distribution of electric charge and 
magnetization as measured by (e, e), the nucleon electromagnetic structure must be folded 
in. This can be achieved using the proton- and neutron electric and magnetic Sachs form 
factors Ge and Gm known from elastic electron-nucleon scattering. In the one-body 
charge- and magnetic form factors, the (appropriate combination of) nucleon form factors 
come in as a multiplicative factor when going from the body form factor (calculated from 
the wave function) to the nuclear electromagnetic form factor. The tacit assumption 
is that the binding of the nucleon in the nucleus does not change its electromagnetic 
structure. One also assumes that the form factors for off-shell nucleons are the same as 
the ones for the free nucleon. For a reliable prediction of nuclear form factors we then 
need to know the nucleon form factors with adequate precision over the range of q of 
interest for nuclear form factors. 

The different nucleon form factors are known from (e, e) experiments to a very different 
level of accuracy: 

The proton magnetic form factor Gmpiq) is well known from many experiments on 
electron-proton scattering. The data reach momentum transfers that go much beyond 
what is needed for a calculation of nuclear form factors (typically 8 fm~^). Several 
parameterizations providing a very good fit to the data are available from Hoehler et al. 
and Mergell et al. 0, ^; other parameterizations can be found in 0, modern references 
are given in p. 

The proton electric form factor Gep{q) is well known at low and medium momentum 
transfers. At the higher momentum transfers, the data up to recently were somewhat 
contradictory. At high q, the magnetic contribution to the e-p cross section increasingly 
dominates, and a determination of Gep{q) via the usual Rosenbluth separation became 
more and more difficult with increasing q. Small systematic errors in the cross sections, 
which depend on a weighted sum of Gl^i^q) and G'^p(g), started to play a big role. 

The situation for Gep has recently improved with experiments that exploit polarization 
observables; when measuring the polarization of the recoil proton in scattering of polarized 
electrons from unpolarized hydrogen, the ratio G ep{q) I G mp{q) becomes accessible 0, 
ni. This ratio is now well measured up to g ~ 9fm~^. This covers the g-range of interest 



for nuclear form factors. The ratio jjLpGep/ Gmp does deviate considerably from the naive 
expectation of 1 as given by the dipole formula. As it turns out, however, the new data 
on Gep{q) largely agree with the parameterization of Hoehler |P which represents an 
often used standard for e-p scattering. This is shown in fig. |l| which displays the recent 
data together with the Hoehler parameterization. Certainly up to g = 8 fm~^, where 
reasonably accurate nuclear charge form factors are available, the proton charge form 
factor Gep, even when used in terms of the venerable Hoehler parameterization, does not 
introduce a significant uncertainty. 

The neutron magnetic form factor Gmn{q), which is important when calculating nu- 
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Figure 1: Ratio of fJ^pGep/Gmp 
parameterization (solid line, 



I [Ty, [TT[] together with the standard Hoehler 
and the more recent fit of Mergell (dotted 



clear magnetic form factors, also is a quantity that is rather difficult to determine exper- 
imentally. Due to the absence of a free neutron target, Gmn has to be extracted from 
experiments on quasi-elastic electron scattering on the deuteron (or more complicated 
nuclei). In the case of inclusive (e, e') quasi-elastic scattering, the subtraction of the dom- 
inant proton contribution, and the corrections for deviations of the reaction mechanism 
from PWIA, lead to large systematic uncertainties. As a consequence, the data showed a 



large scatter |12 



Clean measurements, which do not depend unduly on the theory needed to interpret 
quasi-elastic scattering and to eliminate the dominating contribution from the proton, 
involve coincidence experiments of the type d{e, e'n) where scattering from the neutron 
in the deuteron can unambiguously be identified (save for corrections due to non-IA 
processes). Such experiments are today feasible with the modern high-duty cycle electron 
beam facilities, and we in the following will, as far as possible, only consider this type of 
data. 

Several such experiments have been carried out recently [|T^, ^ at NIKHEF and 



MAMI. The determination of the absolute efficiency of the detector used to identify the 
recoil neutron was measured at the proton-beam facility PSI which allowed one to produce 
intense high-energy tagged neutron beam of variable energy. At very low q, polarization 
observables can also be exploited as the theoretical understanding of the ^He nucleus and 



the reaction mechanism has progressed to the point where polarized ^He can serve as 
a polarized neutron target. In this case Gmn can be obtained from the asymmetry of 
inclusive, single- arm quasielastic scattering W^, 




Figure 2: Data on the neutron magnetic form factor Gmn divided by the 
dipole expression Gd = (1 + g^/0.710(Gey/c)^)~^. The curve represents a fit 
with a continued-fraction parameterization [p!5| . 



The status of our experimental knowledge of Gmn from coincidence measurements was 
clouded by two recent (e, e'n) experiments that claimed to provide accurate data 0, [l^ 
which, however, greatly differed from the ones shown in fig. |^. Both experiments used 
(i(e, e'n) to measure Gmn- Both relied, however, on a reaction with a three-hodj final state 
{d{e,p)e'n or p{e,n)e'n) to produce the "tagged neutron" which was needed to calibrate 
the neutron detector efficiency. Use of a three-hodj final state for tagging — without 
observation of the third particle, the scattered electron — is not possible as a matter of 
principle, and has been quantitatively shown to produce erroneous results |2y, The 
respective results should be ignored. 

Fig. ^ shows the status of today's data on Gmn [0- Up to momentum transfers 
of g ~ 5/m~^ the neutron magnetic form factor is known with excellent precision from 
coincidence experiments [pTS] , pT| , p!5[1 , and these results have been confirmed by experiments 
exploiting polarization observables [17, 16]. At the larger transfers one still has to rely 
on single-arm measurements At large q the exploitation of (e, e') actually gets a bit 

contribution becomes small, so that the (e, e') cross section 

and Gmn- 

On the whole, Gmn is now known with adequate precision over the entire g-range needed 



easier than at low q as the Gep 
then depends mainly on + with comparable contributions from Gmp 
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for the calculation of nuclear magnetic form factors. 



Jourdan et al. |]T5| have parameterized the modern data for Gmn using a continued 
fraction expansion 

G^„ = l/(l + g2.ci/(l + g2.c2/(l + ... 

with the constants Ci,... = 3.260,-0.272,0.012,-2.5,2.5 and q in GeV/c. This fit also 
for the first time allows one to determine to better than 10% the magnetic radius of the 
neutron. It amounts to 0.873 ± 0.011/m (systematic errors included). 

The neutron electric form factor Gen is the quantity that is the poorest known among 
the nucleon structure functions. Only the slope dG eni.^f') / d{q^) of Gen at g = is well 
known from neutron-electron scattering |^3[. The lack of knowledge on Gen can occasion- 
ally really hamper the understanding of nuclear form factors. It does not always do it, as 
the smallness of Gen, which makes its experimental determination so difficult, at the same 
time leads often to a small contribution of the neutron to nuclear form factors. However, 
in interference response functions Gen can play a much more important role. 

The importance of Gen for nuclear form factors can be appreciated by considering 
the charge form factor for an N = Z nucleus. Nucleon finite size leads basically to a 
factor {Genii) + Gep{q)) multiplying the body form factor calculated from the ground 
state wave function. At g = Afm~^ the ratio Gen/ Gep amounts to 0.22 already 0, at 
larger q it could be much larger if Gen{,q) continues to fall as slowly as present data for 
q < 4fm~^ seem to imply. Parameterizations such as the one of Gari and Kriimpelmann 
have predicted large ratios at large q, e.g. Gen/Gep=l.4: at g = 9/m~^ Experiment 



p2| only provides an upper limit of ~ 0.8 (see fig. 0). The uncertainty in Gen is of course 
particularly detrimental for where the neutron charge form factor comes in with 

the double weight. 

The determination of Gen is difficult as there is no free neutron target; in addition, 
the measured electron-neutron cross section is largely dominated by the contribution from 
the magnetic form factor G^n- In the past. Gen therefore was mainly determined from 
electron-deuteron elastic charge scattering which eliminates the magnetic contribution. 
In this approach, the removal of the deuteron structure (deuteron body form factor), the 
subtraction of the proton contribution and non-IA contributions to the reaction mecha- 
nism lead to large systematic uncertainties, however. The most reliable determinations 



were probably those of Platchkov et al. pB] and Galster et al. |2^ which, using various 



calculations of the deuteron wave function and non-IA contributions, provided values of 
Gen up to 4 fm~^. The resulting parameterizations have considerable systematic uncer- 
tainties due to both the deuteron input and the different theoretical results for the non-IA 



contributions |2|, |2|, |0|, |T|. 

The neutron electric form factor can, in principle, also be determined via quasi- 
elastic electron-deuteron scattering, subtraction of the proton contribution and longitu- 
dinal/transverse (L/T) separation. The resulting uncertainties are substantial, however; 
at low q the Gen contribution is comparable to the experimental systematic uncertainties, 
at high q (where it seems that the neutron and proton contributions are more similar in 
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size) the overlap of quasi-elastic and Delta peaks makes difficulties. In fig. ^ we only show 
the results of Lung et al. at very large q. 
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Figure 3: Gen from double-polarization experiments of type d{e,e'n), 
^He{e,e'n) or d{e,e'n) (for references see text). The parameterizations by 
Galster et al. and Platchkov et al. (for the deuteron A{q) calculated with Paris 
and AV14 potentials, respectively) are given by dash-dotted and dashed lines. 
The extrapolation of the parameterizations to large q are shown as dotted 
lines. The Hoehler Gep is also indicated as a dotted line. Also shown are some 
upper limits of G^n from d{e, e') |Q at large q. 



Today, Gen can be determined via polarization observables in (e, e'n) experiments. 
The detection of the knocked out neutron — feasible at CW electron beam facilities 
— largely removes the contribution of the proton. The use of polarized electrons and 
polarized target (or, alternatively, measurement of the polarization of the recoil neutron) 
allows one to measure an interference term G^n ■ Gmn which, once Gmn is known, provides 
a measurement of G^n- 

Given the experimental difficulties of polarized targets or recoil polarimeters, only a 



few results from d{e, e'n), ^He{e, e'n) and d{e, e'n) are available ^21- |^ and the data still 



have fairly large uncertainties. Fortunately, however, these uncertainties are mainly sta- 
tistical in nature, and thus can be improved upon with the progress in polarized target and 
polarimeter technology. The presently available set of Gen data from double-polarization 
experiments is displayed in fig. ^. Clearly, the extent of the data still is far less than what 
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is desirable as an input to the calculation of nuclear form factors of light nuclei. 



2.2 N-N potentials 

An important ingredient to the calculations of the few-body wave functions is the nucleon- 
nucleon (N-N) potential employed. As compared to other interactions occurring in physics, 
the N-N interaction has an unusually rich and complicated structure. It depends on the 
spin and isospin of the nucleons and their relative momentum. Only the long-range part, 
dominated by one-pion exchange, can easily be calculated. The construction of N-N 
potentials is a field that has a long tradition and covers many aspects which we here 
cannot possibly do justice to; so we refer the reader to the original literature (for a review 
see e.g. papers by Machleidt 0). Here, we only can mention a few aspects that are 



relevant for the specific few-body calculations discussed. 

The various N-N potentials are all fit to the N-N (p-p and p-n) scattering data and 
the deuteron binding energy. Typically, the scattering data up to energies of ~350Me\^ 
are employed. At energies higher than 350Me^ inelasticities due to vr-production etc. 
lead to a much more involved description of N-N scattering. While earlier potentials 



43| , |4^ , ^ were generally fitted to the full set of data available, more recent fits 



47| , ^ often have used a restricted data set |Q (see below). 



Most often, the nuclear force is parameterized in terms of a sum of One Boson Ex- 
changes (OBE). These OBE potentials for the longest range component mostly use the 
local version of the one-pion exchange potential. For the intermediate and short range the 
potentials are parameterized in different ways. The Argonne V18 potential developed by 
Wiringa et al. [^] has a one-pion exchange piece for the longest range, but uses essentially 



phenomenological parameterizations for the shorter ranges. The Nijmegen potential of 
Stoks et al. |Q uses empirical masses of existing mesons and meson distributions. All 



potentials are regularized at short distances by e.g. exponential cut-off factors. These 
cut-off form factors simulate the effect of the exchange of heavier mesons which have not 
been taken into account. Modern versions of these OBE potentials are 'charge dependent', 
i.e. they provide different parameterizations for p-p, n-p and n-n. 

The CD-Bonn OBE potential differs mainly in the use of a non-local one-pion exchange 
term, which is derived from a relativistic Lagrangian for meson-nucleon coupling. This 
leads to a non-local force with a tensor component which is substantially weaker off-shell. 
As a consequence, this potential predicts a deuteron D-state probability which is ~0.8% 
smaller than the one predicted by other potentials (the other deuteron static properties 
are very similar). Also, the A=3 binding energy for the Bonn potential is 0.4-0.5MeV^ 
larger. 



Several potentials include more than OBE. The Paris |51| and full-Bonn |^2[ potentials 
include 27r-exchange (or phenomenological fits to it), and thus largely can do away with the 
fictitious cr-meson exchange. The two-pion contribution is calculated using time-ordered 
perturbation theory, or via a dispersion-relation analysis of tt — scattering. Terms like 
the vrp and ttuo exchange have been included as well |^ . The full-Bonn potential is more 



difficult to use in nuclear structure calculations due to its energy dependence, and can not 
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consistently be used for cases of >2 nucleons; as an alternative, the energy-independent 



(but still nonlocal) Bonn-B potential has been derived by Machleidt et al. [41 



The modern OBE potentials fitted to the cleansed N-N scattering data base (see 
below), the V18, CD-Bonn and Nijm-II potentials, achieve a similarly good fit. The 
major difference between them shows up in the calculations of nuclear properties where 
off-shell properties of the force do matter. In particular, a force with a weaker off-shell 
tensor component leads to a stronger binding of the A>3 systems. 

In the past, the potentials were usually fit to the full data base on N-N scattering. 
More recently, the Nijmegen group has performed a fit to the data employing different 
energy-dependent potentials for different partial waves. The of individual data points 
was then used as a criterion to keep the data, or remove them from the data set. A 
large fraction (~30%) was removed. The remaining, cleansed, data base can be fit with 
very good x^, which makes it attractive to those groups constructing N-N potentials and 
therefore has been used to fit the modern potentials mentioned above. 

This procedure of removing data is fine as long as the parameterized energy dependence 
used in the data elimination procedure is correct, and as long as enough degrees of freedom 
have been allowed for. That this is the case is by no means obvious. It can be checked 
by comparing the resulting phases to the ones obtained from single-energy phase shift 
analyses of the full data set. If there are significant differences (as judged by the error 
bar resulting from the single-energy analyses) then one must assume that the multi- 
energy analysis used to reject data was not general enough and introduced a bias. When 
looking at the comparison of phases resulting from single- and multi-energy analyses. 



- / 
y 
-( 

Li 



200 

Elab (MeV) 



o Bugg 

□ Henneck 

+ VPI 

K Bonn 

300 



Figure 4: S-D mixing parameter ei resulting from single energy phase-shift 
analyses (data points, |5^, ^), the multi-energy fit used to cleanse the 



data set (dashed, |5y]), the phases from the CD-Bonn potential (dotdash 
and the 1^18 potential (dotted, [0) fit to this set. 



one finds worrisome systematic differences in some phases. This is true in particular 
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for the ei phase (see fig. 0) |^ which affects the strength of the tensor force and is 
responsible for the S-D transitions that play an important role for few-body systems 
(deuteron electrodisintegration, A=3 magnetic form factors). 

When using OBE potentials in calculations of the wave function of light nuclei with 
A>2, one generally finds an underbinding, which amounts to 0.5-1 MeV for A=3 and 
several MeV for A=4. This indicates that it is also necessary to consider the role of 
the three-body force (3BF). Such a 3BF results in particular from the internal degrees 
of freedom of the nucleon which have been frozen out when describing nuclei in terms of 
interacting nucleons only. 

The longest-range 3BF involves the intermediate excitation of the nucleon to the A, 
and two pions exchanged with two other nucleons. This interaction, originally given by 
Fujita and Miyazawa [0], is attractive in light nuclei. When introducing the A degrees of 



freedom explicitly in the wave function, as done e.g. by P. Sauer and collaborators |[58 
one finds that the attractive contribution is significantly reduced by dispersive effects 
(energy dependence of effective N-N interaction). Other sources for a 3BF need to be 
considered as well. 

Often used are variants of the Tucson-Melbourne force of Coon and collaborators 
|]59| , |60| , |6T[| which are based on 27r and np exchange. At short range these forces are 
phenomenologically adjusted to produce e.g. the correct binding energy of the A=3 
systems and nuclear matter density. Alternatively, the Brazil 3BF is sometimes used ||62 |. 

A somewhat different approach is followed by the Argonne/Urbana group p3| , |6^ 
who also uses at long range a two-pion exchange shape, but at short range uses a purely 
phenomenological spin- and isospin-independent repulsive term, with parameters that are 
again adjusted to fit the binding energies of light nuclei. 

For the time being, these three body forces all contain adjustable parameters. The 
calculation of the 3BF from first principles is still too difficult, as this force not only has 
to account for non-nucleonic degrees of freedom integrated out when describing nuclei in 
terms of nucleonic constituents; it to some degree has to account as well for relativistic 
aspects which are suppressed in the standard, non-relativistic, calculations. 



2.3 Wave functions 

In order to predict the various charge- and magnetic form factors as measured by electron 
scattering experiments, accurate wave functions of the nuclear ground states are the basic 
ingredient. These wave functions can be calculated employing a number of different 
techniques, some of which we will address in other chapters of this review. In this section, 
we only want to discuss in somewhat more detail the approaches that are used in the 
"standard model of nuclear physics". Often these approaches are non-relativistic, but for 
the deuteron relativistic calculations are available as well. The non-relativistic approaches 
calculate the nuclear wave function as a solution of the Schrodinger equation for nucleons 
interacting via one of the modern N-N potentials which reproduce the N-N scattering 
data with a good x^- Calculations performed within this framework have reached a high 
state of perfection, and they have been applied to all of the light nuclei. 
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In this picture of the nucleus, one has aheady made significant simphfications. Rather 
than allowing for all sorts of constituents — nucleons, excited nucleons, mesons, ... — 
one has projected onto a purely nucleonic basis. 

For the deuteron the wave function depends only on the n - p relative distance (relative 
momentum). The solution of the Schrodinger equation for the corresponding one-body 
problem is fairly straightforward. Technical complications arise with the momentum- 
dependent N-N interactions such as the full Bonn N-N potential for calculations in coor- 
dinate space. 

For the three-body system, a number of different approaches to compute the wave 
function within the above-mentioned frame work have been developed. Over the years, 
some of these approaches to calculate ground-state wave functions have been perfected to 
a high degree of reliability p3|-p^ . 

In the Faddeev approach the Schrodinger equation is reformulated in terms of three 
coupled equations separated according to the pair of nucleons that have interacted last. 
The sum of the three equations reproduces the Schrodinger equation. While the separation 
may appear to complicate the description, it actually corresponds to a simplification as, 
in the limit of three identical particles, the 3 equations differ only by a permutation of 
indices. The variables typically used (for the case of coordinate space calculations) are 
the relative distance vector x of two nucleons and the relative radius vector y of the third 
nucleon relative to the center of mass of the pair (Jacobi coordinates), or the corresponding 
momenta. 

The spin- and isospin degrees of freedom of the three nucleons have to be treated 
explicitly such to satisfy the Pauli principle. Often the wave function components are 
decomposed according to the relative angular momentum of pairs of particles, or one 
particle relative to the CM of the other two. Truncations of the space in terms of the 
angular momenta allowed for often are necessary (and can be shown to be justified). 

Accurate calculations today typically keep all channels where the interacting N-N pair 
is in partial waves with j < 4, which requires 34 channels. For scattering states at high 
nucleon momenta, more channels are included in accordance with the larger total angular 
momenta of interest. These Faddeev calculations can be done both in coordinate- or 
momentum space. The former has the advantage that the (for many observables not so 
important) Coulomb interaction can be included more easily. 

In realistic calculations of the A=3 ground state wave function, one wants to include 
the effect of a three body force (3BF). This force is more short-ranged than the vr-exchange 
force, and gives typically a 2% contribution to the potential energy. 

For the momentum-space approaches, it has been shown by the group of Glockle to be 
practical also to solve for the lower energies (<100MeV", say) the Schrodinger equation for 



continuum states |^9[, a great asset when one wants to access the rich physics information 
contained in the two- and three-body break-up channels. 

The three-body Faddeev equations actually can be generalized to four nucleons. The 
corresponding Faddeev-Yakubovsky equations |]7D|] still can be solved although the number 
of channels required for an accurate description of the observables grows very quickly with 
the number of nucleons. 
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During the last years, the Correlated Hyperspherical Harmonics (CHH) approach also 
has successfully been used to describe the A=3,4 bound states as well as the two-body con- 
tinuum channels. The hyperspherical approach involves rewriting the Faddeev equations 
in terms of new coordinates which use (for the case of coordinate space representation) 
the hyperspherical radius which corresponds to the quadratic sum of the usual Faddeev 
radial variables x and y, and the angle between x and y. This hyperspherical approach 
has long been used, without too much success. Only with the introduction in the wave 
function of explicit N-N correlation operators, with parameters determined in a variational 
calculation, has this approach become a successful one. These correlation operators ac- 
count for the strong state-dependent N-N correlations induced by the N-N interaction at 
small N-N distances. Explicitly allowing for these correlations, as done e.g. by Kievsky 
et al. greatly accelerates the convergence of the calculation in terms of the number 
of hyperspherical harmonics that need to be included. 

The most recent addition to the arsenal of methods are Variational Monte-Carlo 
(VMC) calculations. The VMC calculations [0, |73[ start with a trial wave func- 
tion with a Jastrow core that includes single-particle orbitals LS-coupled to the desired 
J,M values. Pair and triplet spatial correlations are introduced from the very beginning. 
The Jastrow core is then acted on by products of two-body spin-, isospin- tensor- and 
spin-orbit correlation operators. Three-body correlation operators are included as well 
according to the 3BF used. The correlation operators allowed for thus contain the same 
operators as occur in the N-N interaction employed. The wave function is diagonalized in 
a basis of different Jastrow spatial symmetry components to project out the lowest state. 

The VMC trial wave functions have typically 20-30 parameters. Using the Metropolis 
algorithm |73|, the parameters are determined by minimizing the energy. Such calculations 
today are feasible up to A=10, and produce ground state energies that are already quite 
accurate (~ 1 MeV). 

In the Greens Function Monte Carlo (GFMC) approach used by Carlson et al. and 
Pudliner et al. [jTSl , Q , one starts with a trial wave function (often the VMC wave function) 



and operates on it with the imaginary time operator e^'^'^'"^"^^, where H' is a simplified 
Hamiltonian and Eq is an estimate of the eigen value, with r being the imaginary time. 
The excited-state components of the trial wave function will then be damped out for large 
r, leaving the exact lowest eigen state with the desired quantum numbers. The expectation 
value of H is computed for the sequence of different 'times' to verify the convergence. 
Differences between the simplified H' and the true H are treated perturbatively. 

These GFMC calculations converge to the essentially exact wave function. Difficulties 
due to the " sign problem" , which limited the times r one could go to, recently have largely 
been overcome (see e.g. 



2.4 Two-body terms 

The work on electromagnetic observables of few-body systems has shown that two-body 
terms of the nuclear electromagnetic operators give important contributions to the observ- 
ables. The explanation of the longstanding discrepancy between theory and experiment 
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of the n+p radiative capture cross section |7|] and the deviation of the A=3 magnetic 
moments from theory provided the first concrete evidence. These exchange terms in- 
fluence both the electromagnetic current and charge operators. The theoretical treatment 
of these two-body terms has made great progress during the past 20 years, as documented 



in several review papers |8y, |30 



The discussion below will start from the assumption that the nuclear system is basically 
non-relativistic, and that relativistic effects can be taken care of by the lowest-order 
expansion in v/c (for calculations in a relativistic framework see sect. p.5|) . Calculations 
of light-nucleus form factors in a fully relativistic frame work at present are still limited to 
the A=2 system; for the other nuclei we still have to rely on the non-relativistic approaches 
involving as complete as practical a set of relativistic corrections. 

We begin by discussing the non-relativistic limit only. The electromagnetic current 
operator j in a nuclear many body system, which we address first, is related to the Hamil- 
tonian H of the system via the continuity equation y • ] + i[H, p] = where p is the charge 
operator of the system. In non-relativistic nuclear physics one usually assumes that the 
Hamiltonian can be separated into a sum of single-nucleon kinetic energy terms T and 
two-nucleon interactions V. As a consequence, the electromagnetic current operator j 
contains a sum of both single-nucleon and two-body operators. Two-body, or exchange, 
operators thus have to be present whenever the commutator between the N-N potential 
and the single-nucleon charge operator is non-zero. This is the case in particular for the 
isospin-dependent and velocity-dependent components of the N-N interaction. Such com- 
ponents occur prominently in the N-N interaction, as the long-range one-pion exchange 
interaction is isospin dependent, and thus leads to an important two-body operator. The 
non-locality of the N-N interaction also leads to a non-vanishing commutator, hence an 



exchange current operator needs to be taken into account 





Figure 5: Contributions to the two-body currents which are model- 
independent in the leading non-relativistic order. The exchange particles V 
and 'p' are vr- and p-like objects. 

Further two- or many-body terms of the electromagnetic interaction would have to be 
considered once many-body forces between nucleons are allowed for. While such forces are 
indeed introduced to account for the binding energies of light nuclei (and nuclear matter) 
the corresponding many-body terms in the currents are still poorly understood. 

The isospin dependent N-N interaction is associated with the exchange of charged 
particles. Within the simplest models, it involves the exchange of tt's and p's. The cor- 
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responding interactions and the associated current operators which satisfy current con- 
servation are known. This allows one to derive the two-body operators from the N-N 

S3, E^. The corresponding terms have been classified as 
In terms of a diagrammatic representation, these 



interaction employed [g^, 
"model independent" by Riska [p4 



terms, accounting for the exchange of vr's and p's, are displayed in fig. ^. By construction 
these terms respect gauge invariance. When dealing with realistic N-N interactions, these 
terms are replaced by exchange terms involving particles that are vr-like (pseudoscalar ex- 
change with quantum number 0~) and p-like (vector- meson type exchange with quantum 
numbers 1^), in accordance with the underlying structure of the N-N interaction [^, |86 



The identification of the vr-like and p-like contributions is not necessarily dictated by the 
potential ||85|, so the qualification of "model independent" has to be understood with the 
quotes. The additional two-body terms then are not necessarily constrained by current 
conservation. 

These terms a priori do not account for the finite size of the nucleon. A simple (but 
not necessarily unique) approach to include the finite size replaces the coupling constants 
by vertex form factors. When using vertex form factors known from N-N scattering, 
additional vertex currents have to be introduced in order to respect current conservation. 
Fortunately, the numerical consequence of these form factors for few-body electromagnetic 
form factors is not great, as their effects in the tt- and p-like terms largely cancel. 

Other two-body terms related to other components of the N-N interaction, such a e.g. 
a velocity dependence, can be and have been included ^ , but are typically of lesser 
numerical importance. 

While the above exchange current operators can be derived directly from the N-N 
interaction, and thus can be regarded as largely model independent, there are others that 
are not constrained by the continuity equation. These terms, addressed in the following, 
then are more uncertain. These " model dependent" terms (in the classification of Riska 





Figure 6: Two-body contributions involving the Delta, and coupling to two 
different mesons. 

HS^ ) are associated with the excitation of virtual intermediate nucleon resonances and 
exchange operators involving more than one type of meson. 

The most important model dependent transition-type exchange current operator in- 
volves excitation of an intermediary A33 resonance (see fig. with the exchanged particle 
being tt- or p-like. Processes involving higher nucleon resonances are suppressed by larger 
energy denominators. While some of the relevant coupling constants are known from 
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experiment, others have to be derived via model calculations. For electron scattering at 
large q, the g-dependence of the couplings becomes a source of uncertainty. In calculations 
of the nuclear ground state wave function that allow for both nucleons and delta's as done 
e.g. by Henning et al. this term is, of course, included (in a non-perturbative way) in 
the wave function and does not appear as an exchange current contribution ||89[| . 

The second type of " model dependent" two-body contribution accounts for terms of the 
type shown in fig. ^a. Such terms, which involve e.g. the pTij transverse isoscalar exchange 
current or the ujtt'j current, are more uncertain due to the poorly known coupling constants 
involved. The pn'y term can, depending on the observable, represent an important source 
of uncertainty at large q, the uir'y and other terms are numerically not so significant. 
In most cases, though, the model-independent terms dominate, and these are well under 
control. 

Coming to the two-body charge operators, the situation is less favourable. In the non- 
relativistic limit, there are no pion exchange contributions to the nuclear charge operator. 
The (longest range) pion exchange charge operator may be looked at as the contact term 
of the type shown in fig. §3. This term, however, is already of relativistic order and 
cannot directly be linked to the N-N interaction employed to compute the nuclear wave 
function. The form of this two-body term (for the exchange of both vr's and heavier 
mesons) still is somewhat uncertain. To derive it one has to systematically reduce both 
the interaction and the electromagnetic current operator from a fully relativistic to a non- 
relativistic description. These contact-type contributions to the two-body charge operator 
are of considerable numerical importance when calculating the few-body CO charge form 
factors. 

When deriving the two-body charge operators by reduction from a relativistic descrip- 
tion of the coupled system of meson and nucleon fields, such as e.g. done by Arenhovel 
et al. pi] , one also has to deal with the boost term which, at very high q, has a 
considerable effect. 

The diagrams involving vr- and p-exchange again need to be modified relative to the 
propagation of bare particles by vertex form factors in order to account for the finite size 
of the particles involved. These form factors can be taken as to be the same ones that 
occurred in the current operators, where they were constrained by current conservation 
and the N-N interaction employed. 

As for the two-body currents, there are contributions to the two-body charge operator 
involving vector mesons, such as the pir'j term (see fig. ^). These processes, as for the 
current operators, are of fairly short range and contribute only at very large momentum 
transfer. 

The situation concerning exchange currents is different when considering the intrin- 



sically relativistic approaches used by VanOrden et al. and Hummel and Tjon 0, |9T[| 
which solve the Bethe-Salpeter equation. If care is taken to insure that the truncations 
necessary in practical calculations do maintain current conservation, then (in PS-coupling) 
no two-body terms appear other than the pion-in-flight diagram (which usually gives a 
small contribution) and the ones due to the type of diagrams shown in fig. p. The piT'-f 
diagram represents the dominant term, and introduces considerable uncertainty. Depend- 
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ing on whether a rather soft vertex form factor or a rather hard one (from e.g. Vector 
Dominance) is used, the deuteron A{q) and B{q) at large q are close to the data, or too 
high. 

3 Deuteron 
3.1 Introduction 

The deuteron, the only bound two-nucleon system, is one of the fundamental systems 
of nuclear physics. Accordingly, many studies, both experimental and theoretical, have 
been devoted to it. Of particular interest today is the degree to which the deuteron can 
be understood as a system of two nucleons interacting via the known nucleon-nucleon 
interaction. 

When addressing, more specifically, the electromagnetic properties of the deuteron, 
the corresponding question concerns the ability to predict the three deuteron form factors 
starting from the calculated deuteron wave function and the nucleon form factors known 
from electron- nucleon scattering. At low momentum transfers predictions and data agree 
quite well when accounting for one-body terms only, at the higher momentum transfers, 
two-body contributions are known to be important. Whether quark degrees of freedom 
do need to be allowed for is still a matter of debate. 

The approach traditionally taken for the two-nucleon system has been to derive the 
nucleon-nucleon (N-N) interaction from N-N scattering data, without inclusion of the 
information from the deuteron (other than the binding energy). The data for elastic 
electron-deuteron scattering then are used as a " check" of our understanding of the N-N 
system. 

The study of the electromagnetic form factors of the deuteron is complicated by the 
fact that the deuteron has spin one. As a consequence there are three form factors corre- 
sponding to the multipolarities CO, Ml and C2, needed to fully describe the electromag- 
netic structure. Without data on polarization observables — the measurement of which 
up to recently presented too many difficulties — only a combination of CO and C2 could 
experimentally be determined, a fact that limited greatly our ability to fully exploit the 
information contained in the deuteron structure functions. This limitation fortunately 
now has largely been overcome. 

The deuteron charge form factor is particularly interesting for the understanding of 
the role of exchange currents. To leading order, (non-relativistic) two body currents are 
absent due to their isovector character. The dominant two-body currents in the isoscalar 
charge form factor then are of rclativistic, v'^/c^, order. These currents still are less than 
perfectly understood. Therefore the CO form factor is a good place to make comparisons 
between the non-relativistic approaches — supplemented by v^/c^ relativistic corrections 
— and the theories based on a fully relativistic approach. 

The deuteron electromagnetic form factors most often are studied in order to check our 
understanding of the two-nucleon system. In parallel, however, the deuteron form factors 
are also exploited to get a better handle on the neutron form factors. In the past, much 
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of our knowledge on the neutron charge form factor Genio) came from precision studies 
of the deuteron structure function A{q). Only very recently, experiments involving both 
polarized electrons and polarized target/recoil- nuclei have allowed us to get access to Gen 
via other observables (see sect. pTTp . At large g, however, Gen is still largely unknown, a 
fact that represents a serious handicap for the quantitative understanding of the deuteron 
charge form factors. 

The deuteron is a very loosely bound system (probably the still best example for 
a "halo nucleus"). Accordingly, many of the deuteron properties of relevance to strong- 
interaction physics concern the large-radius behaviour of the deuteron wave function which 
is well under control. With the electromagnetic form factors at the higher momentum 
transfers g, one does have the chance to access the shorter-range properties of the two- 
nucleon system. This makes the study of the deuteron form factors — especially the CO 
form factor — particularly interesting. 




Figure 7: 
and medium 
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At small inter-nucleon distances, the deuteron actually has quite a bit more structure 
than appreciated when looking at the usual plots of the radial wave functions u{r) and 
w{r) of the S- and D-state, respectively. The strong tensor force at small r leads to 
an interesting behaviour of the deuteron density: for spin-projection = ±1 the equi- 
density contours have the shape of a dumbbell, for mz=0 they have the shape of a donut, 
as displayed in fig. |^ from Forest et al. ||92[ . 
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Traditionally, the calculations of the deuteron electromagnetic structure are compared 
to the structure functions A{q) and B{q) and, more recently, to the data on the polar- 
ization observable T20 that are becoming available. In the present review, we mainly 
compare to the CO, C2 (and Ml) form factors rather than A{q) and T20. This is done for 
two reasons. First, the individual CO and C2 form factors are much more sensitive to the 
physics ingredients than A{q). Second, the use of CO and C2 allows us to more directly 
compare to the A=3 (T=0) and A=4 charge form factors. 

We also want to deemphasize the discussion of the observable T20 for another reason: 
this quantity refers to a specific scattering angle (although the dependence on the angle 
is weak as compared to the one on momentum transfer). It thus is better to avoid this 
dependence on angle altogether and concentrate on CO and C2. Moreover, the physical 
significance of T20 is rather convoluted (even if expressed in terms of the form factors Foo 



and Fii, see ref. |^ and below). We therefore rather use the CO and C2 form factors. 



which today can be determined with adequate precision from the experimental data. 
3.2 Electron scattering data 

Elastic electron-deuteron scattering has been investigated in many experiments, and the 
cross section data |Q - [ |116|| today cover a large range of momentum transfers. Some 
of that data obviously is not very precise, other data, mainly of more recent origin, has 
reached accuracies down to the 1% level. In the analysis to be described below, we will 
employ the full set of world data today available. We discuss below in more detail only 
some of the sets. 

The earliest measurements useable in todays determination of deuteron form factors. 



were performed at Stanford by Friedman et al. ||107|| . These measurements actually 



were aimed mainly at the neutron magnetic form factor, which was studied via magnetic 



scattering from the deuteron. The experiment of Rand et al. ||113|| provided the first more 
extensive set of data on the magnetic form factor measured at 180°. The main result 
of this experiment actually concerned threshold electro-disintegration, the process that 
later on | JL17] | was recognized to be one of the showpieces for the contribution of meson 



exchange currents (see sect. |3^ ). 



The experiment of Elias et al. ||106||, performed at the Cambridge electron accelerator. 



provided data up to 6/m~ . In this experiment, the electron beam was scattered from 
a target placed inside the synchrotron. Elastic scattering was identified by detecting the 
recoil deuterons using a quadrupole spectrometer. These data cover an important range 
in q, where few other data were available up to very recently. Some care is advisable in 
using this data: the deuteron spectra, measured with a combination of scintillators, spark 
chambers and range counters, appeared to show an unidentified constant "background" 
of 20 50% under the deuteron peak. This "background" was subtracted by the authors. 
After closer inspection of the deuteron spectra, and with the benefit of hindsight, it 
would appear that there actually was no significant background; rather, the deuteron 
peak (measured over too narrow a momentum region to allow one to identify a "constant 
background" ) was simply a bit wider than expected. 
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The Bonn synchrotron also contributed to the world supply of electron-deuteron data 



104|| . Using two magnetic spectrometers, the reaction products were detected in coinci- 
dence at both forward and backward angles, in order to perform a selfconsistent Rosen- 
bluth separation. (Note that the last entries in table 2 of ||104]| seem to be inverted). 
Data at the highest momentum transfers were measured by Arnold et al. at SLAC 



96, 101) 112| both at forward and backward angles. In both experiments scattered 



electron and recoil deuteron were detected in coincidence. For the case of the backward- 
angle scattering, the reaction products were detected at zero and 180°, respectively. This 
experiment identified for the first time the diffraction minimum occurring in B{q) at 
q ~ 7fm~^. 

Two recent experiments at Jefferson Laboratory by Abbott et al. and Alexa et al. 
provided additional data P3|, M. These experiments could reach very large momentum 



transfers due to the large beam intensities available at JLAB and the large acceptance 
of the spectrometers. While the hall-A experiment was carried out using the two high- 
resolution spectrometers, the hall-C experiment employed for the recoil-detection the 
magnetic channel built for the T20 polarimeter (see below). In this latter experiment the 
spectrometer acceptance and target length were reduced on purpose in order to achieve 
a well known solid angle and spectrometer acceptance. At the time of the writing of 
this review, there still is a significant discrepancy between the two data sets; there is the 
suspicion that this is due to discrepancies in the determination of the beam energy. 

Of particular interest to the precise determination of the deuteron form factors are the 
sets of data which have reached the highest accuracy. At the low and medium momentum 
transfers, these are the data of Galster et al, Berard et al., Simon et al. and Platchkov 



et al. [0, 100 , TTA , 26]. Precision experiments on the deuteron face extra problems as 
compared to experiments involving heavier nuclei, for two reasons: 

• Achieving a reasonably large target thickness involves in general the use of long liquid- 
deuterium targets. The determination of the acceptance of spectrometers for long targets 
requires a special effort. This poses particular problems for modern multi-element spec- 
trometers which have a complicated acceptance function. 

• Due to the low deuteron mass, the energy loss of the electron due to deuteron recoil 
is large. At the larger momentum transfers the deuteron elastic peak then overlaps with 
the quasi-elastic peak of the target windows. This in general requires to move the win- 
dows outside the spectrometer acceptance, or to detect electron and recoil deuteron in 
coincidence. 

The experiment on e-d scattering carried out by Galster et al. was performed at the 
DESY accelerator using the extracted beam. In this experiment the scattered electron 
and recoil deuteron were detected in coincidence in order to identify elastic scattering 
at the incident electron energy of 2.5 GeV. Galster et al. did a very precise comparison 
between e-p and e-d scattering, with the aim to determine from the deuteron cross section 
the small contribution of charge scattering from the neutron. This determination of Gen-, 
performed up to momentum transfers of 4fm~^, actually has withstood quite well the 
test of time. 



The low-g data of Berard et al. ||100|| were taken using cooled H2 and D2 gas targets. 
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in the range of momentum transfer q = 0.2 0.7/m^^. The experiment essentially mea- 
sured ratios of cross sections relative to Hydrogen; to get cross sections for the deuteron, 
normalization to absolute data on the proton is needed (for small corrections, see 

The data of Simon et al. ||114|| covers the range 0.2 -h 2fm~^. The experiment used 
both a low-temperature gas- and a liquid target, for both Hydrogen and Deuterium. The 
Hydrogen data taken with the gas target and a special small-angle spectrometer served 
as the absolute cross section standard, the data with the liquid targets were measured 
relative to that standard. 

The experiment of Platchkov et al. provided absolute data in the range q = 
0.7 -h 4.5/m~^ measured with a liquid Deuterium target. This experiment reached a very 
high accuracy, This became possible as the spectrometer acceptance was restricted 

to the purely geometric one, which is easily measured. The acceptance of the spectrometer 
for finite target length had been studied in detail using solid targets displaceable along the 
beam direction. The absolute efficiency of the detection system in the focal plane could be 
measured precisely by exploiting the redundancy of the various detector elements. Data 
taken with a liquid Hydrogen target were used to confirm the accuracy of the absolute 
cross section measurement. The Saclay group also provided data on magnetic scattering 
up to momentum transfers of 5.3/m~^ p^, measured at a scattering angle of 155°. 



Some of the deuteron data ||93|-| 116 | listed above have been taken by determining 
the overall normalization using electron-proton scattering. The cross sections given in the 
publications were obtained by normalizing the measured ratios using the best information 
on the parameterization of the proton form factors available at the time. Today, we have 
a set of more accurate and more complete proton data, and a comprehensive fit to all the 
data. It is therefore advisable to renormalize the corresponding deuteron data using the 
fit to todays world data on the proton in order to obtain the most accurate absolute 
deuteron data. 

For a detailed investigation of the form factors and their uncertainties, it is very 
important to account for the systematic errors of the data. Their effect in general is much 
larger than the uncertainty due to the statistical errors. In the publications cited, the 
systematic uncertainties often are not adequately discussed, and several times important 
sources of systematic errors other than normalization [e.g. electron beam energy, beam 
halos) are not even mentioned. For the data sets that have a particular weight in the 



determination of the form factors, one can consult the corresponding thesis works [119 



124|| in order to extract a reasonable estimate for the relevant systematic uncertainties. 



During the last years, it has increasingly become possible to measure not only cross 
sections, but also spin observables. The knowledge of these spin observables is imperative 
if one wants to separate the two contributions of the CO and C2 multipolarities to the A{q) 
structure function. The separation of CO and C2 is of particular interest in the region 
of the predicted diffraction zero of the CO form factor, near q = 4fm~^. To separate CO 
and C2 one needs data involving tensor spin observables. 

Two techniques basically are available to measure such spin observables: 

• At storage rings, one can use polarized, internal deuteron gas targets from an atomic 
beam source. The high intensity of the circulating electron beam allows one to 
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achieve acceptable luminosities despite the very low thickness of the gas target. 



• At facilities with external beams, one can use polarimeters to measure the polariza- 
tion of the recoil deuterons. High beam intensities are a prerequisite as the polar- 
ization measurement, which requires a second reaction of the deuteron, involves a 
loss of a few orders of magnitude in count rate. 

A further potential approach, the use of an external polarized target (such as ND^ at 
very low temperature ||125|) has not yet been employed, as the luminosities reachable with 



these targets were not sufficient. This in part is also related to the difficulties of making, 
instead of the usual vector-polarized targets, a tensor-polarized target as needed for the 
separation of CO and C2. 

The pioneering experiment on e-d scattering involving polarization observables was 
performed by Schulze et al. ||126|| at Bates. The polarimeter employed made use of the 



^He{d,pYHe reaction, which represents an efficient analyzer of the tensor polarization for 
deuterons with very low momentum. This experiment provided data near q ~ 2fm~^. 



^H(j2p)n 

Ed= 200 -- 400 MeV 




Figure 8: POLDER deuteron polarimeter used in the experiment of Abbott et 
al. |^3[. The deuteron track is measured before the d + p-^p + p + n reaction 
using multi-wire proportional chambers. The two protons are detected using 
the hodoscope counters. 



The first experiment with an internal electron beam and a gas-jet polarized target was 
performed at the VEPP-storage ring in Novosibirsk ||127|| . The data cover the region near 
lfm~^ momentum transfer. The follow-up experiment at VEPP |[L28|| reached 3/m~^, 
but had rather large statistical uncertainties. 
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These early experiments did not yet cover the region where the separation of CO and 
C2 is expected to provide new information, but they showed that — with a suitable 
increase in luminosity and polarimeter efficiency — the determination of polarization 
observables would be practicable. 

Higher momentum transfers could be reached using external beams and different types 
of polarimeters. The experiment of The et al. |[L29|| , performed at the Bates accelerator, 
used a polarimeter based on d — p elastic scattering. With this setup it became possible 
to reach the region of the zero crossing of the CO form factor at 4.4/m~^. 

Further measurements using an internal tensor-polarized gas target were carried out 



at the AMPS storage ring ||130| , |131|| . The deuterons were produced using an atomic-beam 
source which injected the deuterons into a storage cell traversed by the electron beam. 
Scattered electron and recoil deuteron were measured in coincidence. This experiment 
provided T20- values up to a momentum transfer of 3.2/m~^ 

The highest momentum transfers were reached in the recent experiment performed at 



JLAB by Abbott et al. ||9^. With this experiment the data could be extended beyond the 
diffraction maximum of the CO form factor, reaching 6.7/m~^. This became possible with 
the construction of a polarimeter optimized for the high deuteron momenta, a polarimeter 
based on the d + p —^p + p + n break-up reaction depicted in fig. |[ This experiment also 
resolved a worrying discrepancy between the data and theory in the 4-5/m~^ region where 
the normally most successful calculations seemed to lie significantly above the previous 
data. 
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Figure 9: World data for the tensor analyzing power T20 for electron-deuteron 
scattering as function of momentum transfer (the slight dependence on scat- 
tering angle has not been removed). 



In fig. 1^ we display the T20 data which today are available. Additional data have 
been measured for T21 and T22. These quantities, however, are not very constraining in 
separating the multipolarities CO and C2; rather, they serve as consistency checks. 
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3.3 Experimental form factors 



Here we use the world data on electron-deuteron scattering up to high momentum transfer 
to determine the various form factors. Some 340 data points on e-d scattering are available 
for momentum transfers below 8fm~^ (the higher transfers are ignored for the moment 
as the experimental information is more fragmentary). 

As discussed in more detail in section ^.4| , the experimental cross sections are first 
converted to effective PWIA cross sections by removing the Coulomb distortion effects. 
This is done by multiplying the experimental cross section a{E, 9) with the calculated ratio 
ci.BorniE, 9) / a2.Born{E ^ 9). Thcsc a^^^^ then can be interpreted in terms of PWIA. 

The deuteron cross sections measured in the various experiments [p5|- [|116|| contain 



contributions from the form factors of multipolarity CO (charge monopole), Ml (magnetic 
dipole) and C2 (charge quadrupole), which in turn determine the two structure functions 
A{q) and B{q): 

a{E, = aMottiE, 9) [A{q) + B{q) tg\9/2)] 

where aMott includes the recoil factor, with 

A{q) = F^q) + {MjQf-r'FUq) + i^^f^f^il + r)FUq) 



In the above equations, M are the masses, /i the deuteron magnetic moment in units of 
magnetons, Q the quadrupole moment and q is the four momentum transfer {'JS). All the 
F-form factors are normalized to -F(O) = 1. 

Up to recently the study of the form factors in general was limited to the A{q) and 
B[q) structure functions. With the polarization data that became available during the 
last years ||126| , |127| , p,28| , |116| , |130|| , one also can perform an analysis of the world data 



in terms of the form factors Fcoic[), F]\ji(q), -Fc2(q')- The tensor polarization observables 
give a handle to separate the contributions of CO and C2. In particular 

T20 = -ilrFcoFc2QM^, + Ir^F^^.Q^M^, + r(l + ^(1 + r) tg^-^Fl,{^^^f)l ^2{A + Btg' 

contains the interesting interference term CO ■ C2. The other tensor observables, T21 and 
T22, are not so useful. 

The separation of Ai^q) and Bicj) (or -Fco((?), FMi{q), Fc2{,<i)) in the past often has 
been done by the experimental groups providing the cross section data. In most cases, 
the individual experiments measured cross sections that were dominated by either A{q) 
or -B(g), so that the "Rosenbluth separation" could simply be done by subtracting the 
non-dominating structure function using the best values known at the time. 
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Here, we rather concentrate on a more reliable approach which involves a fit of the 
unseparated world cross sections with a parameterization of A{q) and B{q) simultaneously. 
This allows one to use the full set of data today available to determine all form factors. 
The "Rosenbluth separation" thus is done implicitly during the fit of the data. As the 
parameterizations of A{q) and B{q) provide continuous functions, this approach also takes 
care of the problem that rarely the forward-angle and backward-angle data are available at 
exactly the same q. For the standard Rosenbluth separation this usually implies somewhat 
obscure interpolations or extrapolations. The same reasoning applies when determining 
the three form factors CO, Ml and C2 by adding the polarization observables. 

In order to analyze the deuteron data, the various form factors are parameterized 
using a flexible functional form. From these form factors, one calculates cross sections in 
PWIA, and fits the parameters to the data (already corrected for Coulomb distortion). 

We here discuss the results obtained with the Sum-of-Gaussians (SOG) parameteri- 
zations jl32|| for the F's or, alternatively, A{q) and B{q). The form factors are written 



as 



Qi ( / n N 2i?j^ sin{qRi 



+ 



For the comparison to theory mainly the g-space version is of interest a priori. In co- 
ordinate space, this parameterization for -Fco(q') corresponds to the Fourier transform of 
a density p(r) written as a sum of (symmetrized) gaussians that are placed at arbitrary 
radii with amplitudes Qi that are fitted to the data, and a fixed width 7. 

The choice of this parameterization is governed by the following considerations. To 
extract form factors not biased by the choice of the analytical form, one needs a param- 
eterization that in principle is totally general, with restrictions that can be justified on 
physical grounds. In the limit 7 — 0, i?„ — 00 (hence n — > 00) the SOG is a perfectly 
general basis. For a fit of the data, one introduces two restrictions: 1) The gaussians are 
given a finite width 7. In both Pchij") and Pmir) we do not expect any structure that is 
significantly smaller than allowed for by the proton size. Accordingly, 7 is defined via the 



rms-radius 7 y3/2 = 0.6 fm. It has also been verified that with this value the form fac- 
tors from several theoretical calculations can be accurately represented up to the largest 
q. This finite width essentially imposes that the form factors on average fall no slower 
than the proton form factor. 2) The gaussians are placed at radii Ri < Rmax = 10/m. 
This is justified given the fact that, from independent knowledge on the behavior of wave 
functions at large radii, one can easily specify the radius at which the tails of densities give 
no significant (< 10"'^) contribution to F{q). A finite Rmax imposes that the form factor 
cannot oscillate in g-space more quickly than with a minimal wave length Aq cx 1/Rmax- 

The data have been fitted with the SOG parameterization, and 10 free parameters (for 
r < 4/m) for each structure function (for the treatment of the region 4 < r < 10 fm see 
sect. p.4|) . The data set for the deuteron is thus fitted with a total of 30 free parameters 
(20 for the fit of A(g), B{q)). 

The statistical error of the fit is easily calculated using the error matrix. The system- 
atic errors, which in general are the dominating ones, have been evaluated by changing 
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each individual data set by the quoted error, and refitting the complete data set. The 
dominant systematic errors are the ones due to overall normalization, and knowledge of 
the electron energy. The changes due to systematic errors of the different, independent, 
sets of data are evaluated separately, and added quadratically. The total error then is the 
quadratic sum of statistical and systematic errors. 

With the SOG parametrization both the form factors Fco, Fmi, Fc2 and the structure 
functions A{q), B{q) have been extracted from the data. (Numerical values are available 
upon request). For the former quantities, the separation is limited to a maximum mo- 
mentum transfer of ~ 7fm~^, due to the limited set of T20 data available. 

The of the fit is quite satisfactory, given the very diverse origin of the data. When 
fitting the three form factors to the data where the uncertainties include the statistical 
errors only, the amounts to ~610 for 373 data points Q (In this fit the large-r behaviour 



was constrained as described in sect. p.4|) . When adding quadratically the systematic 
error, as is often done, the of this fit amounts to 485. A large contribution to the 
(~110) comes from the data of ref . ||106|| which, as discussed in sect. |3.2| , suffer from a 



probably incorrect background subtraction. 

The form factors resulting from the fit, together with their total error (statistical and 



systematic) will be used in sect. p75| when comparing to the theoretical calculations. 

Traditionally, the deuteron electromagnetic structure has been discussed in terms of 
the A{q),B{q) structure functions or the Fco, Fmi, Fc2 form factors. As pointed out by 
Forest et al. [Q, it may be more elucidating to employ, for the charge form factors, the 
quantities which (in lA) correspond to the Fourier transforms of the deuteron density 
in welldefined states of the projection = 0, ±1 of the deuteron spin, Fqo, Fu. (To 
distinguish these form factors from the usual ones, we use subscripts to indicate the initial 
and final spin projection 0, ±1). These form factors are linear combinations of the usual 
CO, C2 ones (the corresponding magnetic form factor Fqi, which describes the transition 
from the = to the niz = ±1 state is equivalent to Fmi): 

Fco = 1/3 Foo + 2/3 Fu Fc2 = 2 Foo/Qq' - 2 F^/Qq^ 

We give in fig. ^ plots for these quantities as the uncertainties of the Fqo, Fu are not 
easily computed from the ones of Fco, Fc2 shown in other figures of this review. 

At large q, the form factor Fu gets its largest contribution from the two peaks of pi(r) 
occurring at distance r = d from the deuteron center of mass CM (see fig. 0). If these 
peaks would have very narrow width, the form factor would have a cos(gc/)-dependence, 
with the first sign change occurring at qd = vr. The position of the first zero of Fu thus 
gives information related to the radius r where the maximum density occurs. 

When approximating the deuteron density in the = state as a disc of thickness 
t, the corresponding form factor would have a sin{qt/2)/{qt/2) dependence. The first 
diffraction zero thus gives information on the thickness t; the presently available form 
factors do not yet allow to determine it directly. 

"'^Due to the presently still unresolved discrepancy between the data sets of [BSl M mentioned in 



sect. 3.2 we have included only the set 
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Figure 10: Charge form factors of the deuteron for spin-projection = 
and rriz = 1. The error bars include statistical and systematic uncertainties. 
The experiment is compared to selected theoretical calculations (for discussion 
of the calculations see sect. |33| ). 



When ignoring the (rather small) magnetic contribution to T20, this observable is in a 
fairly simple way related to the = 0, ±1 form factors: 

The minimum of T20 thus occurs when Fii=0, while for the maxima Foo=0. The minima 
and maxima thus occur at those values of q where the recoiling deuterons are only in the 
ruz = and = ±1 states, respectively. 

Similar considerations apply to the magnetic form factor |Q. The zero of the mag- 
netic form factor, which experimentally is located near 7fm~^, gives information on the 
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thickness of the torus, which amounts to about 0.9/m. 



In fig. |T0|, we give the experimental results from our analysis of the world data and 
compare to some theoretical form factors for the cases = 0, ± 1; a more detailed 
discussion of the theoretical calculations is given in section pTS] . 

In the above figures, as in most of the other figures in this review, we plot the form 
factors as a function of q, and not, as is often done, g^. The use of leads to an 
overemphasis on the highest g's where the uncertainties of the data generally are largest. 
To the degree that two-body effects can be ignored, the densities corresponding to the 
form factors are Fourier transforms involving g as a variable, and not g^. 

3.4 Radius 

The electromagnetic properties of the deuteron at low momentum transfer could be hoped 
to be accurately predicted as non-nucleonic degrees of freedom or relativistic effects — 
aspects that are not yet under perfect control — should be unimportant. The form factors 
at very low g are dominated by the parts of the deuteron wave function where the two 
nucleons are far apart, and the properties of the deuteron should be determined by the 
known N-N interaction and the known nucleon form factors. 

For these reasons, the deuteron rms-radius has been a favourite observable to compare 
experiment and calculation. The theoretical calculations of the rms-radius are particu- 
larly reliable as the calculation is largely independent of the particular nucleon-nucleon 
potential used; for a very broad class of nucleon-nucleon potentials the radius depends 
essentially on the binding energy and the well known n-p scattering lengths or, alterna- 
tively, the known asymptotic norm. A comparison between calculation and experiment 
therefore promises to be particularly constraining. 

Over the years, many theoretical calculations aiming at the deuteron rms-radius have 
been performed (for a review see ||133|| ). Much of the insight from theory, and an analysis 
of the experimental data, have been put together by Klarsfeld et al. ||134|| , who discovered 



a disturbing discrepancy: the rms-radius derived using nucleon-nucleon potentials was 
0.019±0.003 fm higher than the one obtained from electron scattering data. The de- 
tailed analysis of Wong ||133|| confirmed this finding. Given the failure to find a plausible 
mechanism to explain this discrepancy, this led to a rather confused situation. 

This discrepancy has triggered several authors to look more closely into potential 
corrections. In particular, Buchmann et al. and Herrmann et al. have studied the effects 
of meson exchange currents |^ and dispersive corrections ||135|| in much greater detail 



than had been done previously. The effects found, however, were quite minor in terms 
of a change of the rms-radius, and go in the direction of increasing the discrepancy. The 
modification of the deuteron wave function due to other aspects that are not well under 
control, such as the energy- dependence of the N-N interaction off-shell |r^, 137 |, have 
been studied, also with little success in explaining the discrepancy. 



It was finally found | |138|| that much of the discrepancy originated from the fact that 



the deuteron data always were analyzed in Plane- Wave Born Approximation (PWBA), i.e. 
by neglecting the Coulomb distortion. Although Coulomb distortion is indeed small as 
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Za ~ 0.01, the distortion effects are significant at the level of precision the comparison of 
radii from various sources has reached today. They should be taken into account not only 
for the deuteron. Also for the proton it turns out that, when considering the rms-radius. 



Coulomb distortion makes a non- negligible difference, as pointed out by Rosenfelder ||139 |. 

In order to calculate the Coulomb distortion, ref. [ |138|| used the second order Born 
approximation. The series in Za is expected to be very accurate for the Za ~0.01 of 
interest for Hydrogen and, as shown in ref. |135|, there indeed is no significant difference 



between this approach and the exact partial-wave analysis cross section. 

The cross section in second order Born approximation [|140|] can be written, for scat- 
tering angles not very close to ^ = 180°, as 

02.Born ^ CTMott 1^(1 + R) 

where F{(i) is the elastic form factor, q = ki — kf is the momentum transfer, fcj is the 
incident electron momentum and kf is the final electron momentum. R is the Coulomb 
correction factor and is given by a principal value integral as: 
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where k = Iki 



\kf\ and where we have used the abbreviation A = q/p = tan with p = 
ki + kf. This integral, which in its original form is rather difficult to calculate numerically. 



has been simplified ||138|| by regularizing the principal value integral by subtracting the 
singular part of the integral. The resulting integrals are numerically well behaved and 
can be calculated easily for any form factor given in analytic or numeric form. This 
calculation | p,18| ] is used throughout this paper to convert the experimental data to PWIA 
form factors. 

At the low q of main interest for a radius-determination, the dominant effect of 
Coulomb distortion actually is not the one due to the familiar difference q ^ qeff (which 
has little effect indeed). Rather, it is the change in cross section given, for a point nucleus, 
by e.g. McKinley [ [141|| . For a point nucleus, the Mott cross section then reads 



O'Mott 



Ze^y cos\9/2) 
2E J sin\e/2) 



1 + 



TiZ sin{e / 2){l - sin{d/ 2)) 



137 



cos^{e/2) 



The additive term proportional to Z/137 is the one that is mainly responsible for the 
Coulomb distortion at the low g's as it influences the finite size effect which is given by 
the difference between the experimental cross section and the point nucleus value. 

The deuteron is special in the sense that the density extends to rather large radii, given 
the low binding energy. The long tail has caused considerable difficulties in the past, as 
it influences the deuteron form factor at extremely low g, where accurate values for the 
difference of the form factor to the point nucleus value 1 are hard to measure. Fortunately, 
at these large distances, the shape of the deuteron density is easily calculable. Outside 
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the range on the N-N force, the wave functions u{r) and w{r) have an analytic form that 
is well known (see e.g. ||142|| ) that depends only on the deuteron binding energy. For the 
determination of the deuteron radius, this shape for radii r > 4/m, where the distance r 
refers to the distance of the nucleons to the deuteron center of mass, can be imposed. 

The fit of the world data for q < 8fm~^ has been performed as described in section |3.3| . 
For the rms-radius of the deuteron, ref. ||138|| finds 2.130 fm, with a random uncertainty 
of 0.003 fm and a systematic uncertainty of 0.009 fm. The results of the fit is reported 
in table |l|, together with results from optical isotope shifts and deuteron wave functions 
to be discussed below. 

The comparison to other results is made on a basis of the quantity that comes closest 
to the charge rms-radius of the non-relativistic two-nucleon system. We remove from 
the quantity directly measured in (e, e) the effect related to non-nucleonic degrees of 
freedom (two body currents) as far as presently possible. We also remove the effects due 
to the (virtual) excitation of internal degrees of freedom of the deuteron by the electron 
(dispersion corrections). In the same spirit, we remove from the charge radius measured in 
optical transitions the contributions of meson exchange currents and nuclear polarization. 
To compare to the nuclear size calculated by solving the Schrodinger equation for given 
N-N potentials, we add to the "matter rms-radius" (the expectation value of of the 
wave function) the contribution of the proton and neutron charge radii, and the Darwin- 
Foldy term. This makes the various radii comparable; the radius we discuss is closest to 
what could be considered as "charge radius" in (e, e) in Impulse Approximation. 

The dispersive effects — corresponding to a two-step scattering process with excitation 
of the deuteron in the intermediate state — have recently been studied by Herrmann and 
Rosenfelder | 135 | who take into account the Coulomb excitation only and use an S-wave 
separable potential (Yamaguchi) to calculate the deuteron wave function. When analyzing 
the data of Simon et al. ||114|| , they find a change of the rms-radius of -0.003 fm when 



correcting the data for the dispersive effects. This calculation gives a significantly smaller 



effect than a previous estimate |p.43|| , but is much more reliable. The contribution of 
non-nucleonic degrees of freedom have been studied in great detail by Buchmann et al. 
When correcting the experimental data for the two-body effects, these authors find 



a change of the rms-radius of -0.001/m, with a fluctuation of .001/m depending on the 
approach used. The estimate for the contribution of 6-quark components is much smaller, 
and presumably also much more uncertain. 

For the optical isotope shift, the accuracy has recently greatly improved, as a conse- 
quence of progress in the area of two-photon spectroscopy on hydrogen and deuterium. 
The atomic transition energies are now known with much higher precision, and a num- 
ber of additional higher-order QCD corrections have been calculated by Pachucki et al. 
144| , [L45|| . We here quote the analysis of these results by Friar et al. ||137|| . Using for 



the proton charge rms-radius the standard value of 0.862/m correcting for the small 
nuclear polarization effect and using, for consistency, for the two-body effects, the value 
as calculated in we find 2.1316 ± 0.001 fm for the deuteron rms-charge radius. 

The radius of the deuteron as determined from N-N scattering has remained rela- 
tively stable since the work of Klarsfeld et al. [|134|| . Friar et al. have also updated the 
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determination of the deuteron radius starting from the nucleon-nucleon potentials and the 
known asymptotic normahzation of the S-state wave function. This yields, when adding 
the remaining electromagnetic contributions due to proton and neutron intrinsic charge 
distribution and the Darwin- Foldy term, 2.1286 ± 0.002/m. With a small further cor- 
rection resulting from the usually neglected proton-neutron mass difference, we get the 
value listed in table |l|. 



Source 



rms-radius {fm) 



reference 



fe, e) world data 



2.130 ±0.010 



(e, e) -disp. -MEC 
isotope shift -pol. -MEC 
N-N scattering data 



2.123 ±0.010 
2.1316 ±0.001 
2.130 ±0.002 



i TTsj |T35| , H 

29|, pTj PI 

m 



Table 1: Charge rms-radii from different sources. Corrections 
for dispersive effects (disp), two-body currents (MEC) and nuclear 
polarization (pol) have been applied to make the three last entries 
comparable. 



The comparison of the second and fourth entry of table |I| shows that there is good 
agreement between the radius determined from (e, e) and N-N scattering. The radius 
coming from the optical isotope shifts (third entry) is within errors compatible with the 
one from electron scattering, but lies slightly above the value deduced from N-N potentials 
and the deuteron asymptotic norm. Overall, we find quite a consistent picture, as the 
three sources on the deuteron size give results which are quite close. 

Occasionally, the equivalent of the charge rms-radius for the Ml part is of interest. 
This quantity is not really a "radius", but it does describe the g^-dependent term of the 
Ml form factor at g = 0. This "magnetic radius" amounts to 2.072 ±0.018 fm. 

3.5 Comparison to theory 

The theoretical understanding of the deuteron form factors involves a number of issues. 
We first list the main questions before discussing in more detail selected calculations and 
their comparison to the experimental results. 

The impulse approximation (lA) calculation of the form factors in a non-relativistic 
frame work is straightforward as the Schrodinger equation can easily be solved for a given 
nucleon-nucleon (N-N) potential. Uncertainties originate from the N-N potential employed 
(mainly its off-shell properties related to the non-locality of the potential) and the nucleon 
form factors used (mainly the neutron electric form factor Gen)- For the deuteron CO form 
factor, the sum of neutron and proton Gen and Gep comes in as a multiplicative factor; 
uncertainties in Gen then directly propagate in the charge form factor. 
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It may be instructive to consider the lA-relation between the S- and D-state radial 
wave functions u{r) and w{r) and the form factors. The form factors are given by the 
integrals 



POD |- -, 

Fco{q) ^ 2Gl{q) u{rf + w{rf jo{qr/r) dr 

Jo L -I 



Fc2{q) = rw{r) [u{r) -w{r) /2V2\ Uqr/2) dr 

Q ^ 



q'Q 



FmM = - [Gl{q)Fi{q) + 2G^(g)F,(g)] 
A* 

where the Ml form factor is split up into the contributions from the intrinsic magnetiza- 
tion and the spin-orbit term 



3 /'°° 

Fi{q) = 7; / ^irf Uo{qr/2)+j2{qr/2)] dr 
2 Jo 



F,{q) = (\u{rf - w{rf/2\ jo{qr/2) + w{r)/V2 [u{r) + w{r)/V2\ ^2(^^/2)) dr , 
where the isoscalar nucleon form factors are given by 



Gl{q) = [GeM+Gen{q)]/2, 



G'miq) 



[Gmp{q) + Gmn{q)]/2. 



and r here refers to the n — p distance. The D-state contribution affects, besides the C2 
form factor, also significantly the Ml form factor due to the S-D transition, by shifting 
the lA diffraction minimum by about 2fm'^ to larger q and decreases the height of the 
diffraction maximum. The contribution of w in Fco is a modest increase of the height of 
the diffraction maximum, the effect of Fi in Fmi is small. 

In order to quantitatively understand the deuteron electromagnetic structure, the 
contribution of two-body terms needs to be accounted for. In order to maintain cur- 
rent conservation, the two-body currents should be consistent with the N-N interaction 
employed. The calculation of such two-body terms has become possible during the last 
years. Uncertainties arise mainly from those two-body terms (such as the 7rp7-diagram) 
that are not constrained by current conservation. The deuteron is special in the sense that 
the leading-order nonrelativistic MEC from vr-exchange are absent due to their isovector 
character; thus relativistic contributions are expected to be more important. 

In a relativistic frame work — which is needed if one wants to understand the deuteron 
electromagnetic structure up to the largest momentum transfers — the calculation of the 
form factors gets more involved and different approximations arc used by different authors. 

Below, we compare the experimental results to a number of theoretical calculations. 
We first address calculations that assume that the basic dynamical content of the deuteron 
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Figure 11: Comparison of experimental CO form factor and calculations of 
Henning et al. (top), and vanOrden et ai, Carlson et ai, Hummel and Tjon 
and Arenhovel et al. (bottom). For references see text. 



is non-relativistic in nature, and that relativistic effects can be described as corrections in 



a f /c-expansion ||146|| . In these calculations, the N-N interaction is one of the standard po- 
tentials with parameters fitted to the N-N scattering data. The two-body terms are then 
constructed to be, as far as possible, consistent with the meson exchanges of the poten- 
tial. We subsequently address calculations which start from the Bethe-Salpeter equation 
representing a summation of multiple-scattering series written in terms of Feynman dia- 
grams. This approach is covariant, but the solution of the Bethe-Salpeter equation usually 
requires truncations which can lead to problems in maintaining Lorentz covariance and 
current conservation. 

A large number of calculations are available | 



147, 148, 149, 150, 151, 152, 47 



so 
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Figure 12: Comparison of experimental Ml form factor and calculations of 
Henning et al. (top), and vanOrden et ai, Carlson et ai, Hummel and Tjon 
and Arenhovel et al. (bottom). For references see text. 



we here discuss only a few representative ones which have appeared more recently. 

The Hannover group |153| , [154| , |155| , p.56|| has performed calculations for various 
N-N potentials (Paris, Bonn-A,B,C,..). For the nucleonic form factors the authors use 
the Hoehler parameterization. For the calculation of the two-body terms, they use the 
approach developed by Adam et al. who employ a modified S-matrix method to derive 
the two-body contribution including the leading order relativistic corrections. The elec- 
tromagnetic nuclear currents are derived from the exchange of all mesons which the Bonn 
potentials involve. Thus, the two-body contact-, pair- and retardation terms are consis- 
tent with the underlying N-N potential. Boost corrections in the Ml are not included. 
The p7r7 term is added separately. 
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The results of these calculations are reported in the upper parts of figs. In 
order to illustrate the two-body contribution, the lA result is also shown for one particular 
N-N interaction. The comparison to the full calculation shows that for the CO and Ml 
form factors the contribution of two-body terms is very important. 



The calculations of Carlson and Schiavilla are based on the Argonne V18 N-N 
potential. The one-body operator uses the standard Hoehler form factors. The two- 
body current has '"model-independent" pieces, related to the vr-like and p-like exchange 
terms in the N-N interaction; these are largely derived from the N-N interaction and 
are constructed to satisfy current conservation. The mo del- dependent part of the two- 
body current is mainly associated with the A-isobar. By means of transition correlation 
operators AA (and, for A > 2, NA) components in the wave function are created. The 
contribution of the piT'y and ujtt'j exchange charge operators is more ambiguous as it 
requires the introduction of phenomenological cut-off form factors. The results of this 
calculation are also reported in figs. |TT]-|T3[ 



In the figures we also show the results obtained recently by the Mainz group |31[]. This 
calculation starts from a system of coupled nucleon and meson fields, and, by means of 
the Foldy-Wouthuysen transformation, derives the non-relativistic limit including all the 
leading order relativistic contributions. For the exchange currents the terms consistent 
with the Bonn OBEPQ-B potential, which is used to calculate the wave function, are 
used. The lowest order piT'y contributions are also included when calculating the Ml form 
factor. The dipole nucleon form factors (which do not well reproduce Gep at the larger 
g's, see above) together with the Galster Gen are employed. 

This calculation finds, for the CO form factor, a sizeable effect of the boost corrections 
(not included in ^^), which are however largely canceled by the vr-two body term. The 
boost contribution for the Ml multipolarity, together with the vr-two body-term, leads to 
a significantly too large Ml form factor at the larger g's. 

The calculated two-body contributions obviously go a long way in closing the gap 
between the lA results and experiment. Some differences do remain, however. This is 
also seen in the g = observables. For instance, the lA deuteron magnetic moment for the 
V18 potential is 0.847 magnetons, with the two-body terms added Wiringa et al. find 
0.871, while experimentally the magnetic moment is 0.857. Similarly, for the quadrupole 
moment, the lA value of 0.270 fm? is increased by two-body contributions to 0.275 only. 



while experiment finds 0.286. For a discussion of the rms-radius, see sect. p.4[ Other 
types of calculations find similar differences, the origins of which are not understood. 

We next address calculations that belong to the second class of approaches mentioned 
at the beginning of this section. Relativity is imposed from the very beginning by requiring 
that currents and interactions are consistent with Poincare invariance. Several groups 



have published results belonging to this category ||157|| - ||161| , [p0|] . We here can discuss 
only some selected ones. 



The calculation of Hummel and Tjon ||91| employs a quasi-potential approximation to 
the Blankenbecler-Sugar equation. The description of the nucleonic and mesonic struc- 
ture is performed within the one-boson-exchange (OBE) model, with the various coupling 
constants fitted to the experimental phase shifts of N-N scattering with the same rela- 
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Figure 13: Comparison of experimental C2 form factor and calculations of 
Henning et al. (top), and vanOrden et ai, Carlson et ai, Hummel and Tjon 
and Arenhovel et al. (bottom) For references see text. 



tivistic frame work. For the results shown here, the Hoehler nucleon form factors have 
been employed. To the form factors resulting from this relativistic lA (RIA) calculation, 
the contributions of the puj and ue'y MEC diagrams have been added, with form factors 
taken from the vector dominance model (VDM). 

In figs. p!l|-p!3| we also display the results of the calculations by vanOrden, Gross and 
Devine ||9^. These authors also start from the Bethe-Salpeter equation, which has been 
reduced to a quasi-potential equation by assuming that one of the nucleons is on mass 
shell. This calculation again is Lorentz covariant and gauge invariant. The parameters of 
the OBE interaction have directly been fit to the N-N scattering data. The calculation 
includes the contribution from the pn'-f exchange current, with vertex form factors coming 
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Figure 14: Comparison of experimental A{q) and B{q) structure functions 
and calculations of Henning et ai, vanOrden et ai, Carlson et al. For references 
see text. 



from model calculations. Contrary to most other calculations cited above, here the dipole 
nucleon form factors have been used, together with a particular prescription for the off- 
shell nucleon form factors. 

It is found that wave function components which occur only in the covariant descrip- 
tion, e.g. the p-state contribution, give a rather small contribution, even at large q. The 
main effect of the fully relativistic description then must be assumed to come from the 
use of a fully relativistic electromagnetic operator, as has been found in refs. [[L62| , |163|| . 



When comparing the results of the different calculations, it becomes clear that the 
C2 form factor is not very sensitive to the dynamical ingredients of the calculations. The 
same is true for A{q) which, over a significant part of the g-range, is dominated by the 
C2 multipole form factor. The CO and the Ml multipolarities are the ones that allow for 
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the most sensitive test of our understanding of deuteron structure and two-body effects. 

When considering all three form factors, the calculations using the Bonn-B, Paris and 
V18-potentials yield the best overall agreement with the data. The deviations of the 
Bonn- A curves (C2 form factor falling too quickly, Ml diffraction zero at too low a q) are 
characteristic for a wave function with too large a " defect" at small r in the D-state wave 
function (leading to a low D-state probability). 

The role played by the neutron electric form factor has been investigated several 
times, most recently by Plessas et al. [Q. As pointed out in sect.[2?T| the deuteron charge 



observables F^q, ^iid A{q) basically contain a factor {Gep{q) +Gen{q)Y- Large values 
of Gen{q) at large g, as given e.g. by the parameterization of Gari and Kriimpelmann 



2J] can lead to a large factor. For this parameterization the ratio of the predicted Ge 



to the experimental Gep at g = (the upper limit of the charge form factors shown 

in figs. |ll],|13]) is Gen/ G^p^ 1.2, while the Hoehler parameterization yields 0.3; the data 
only provide an upper limit of ~0.8. 

In order to judge the agreement between theory and data at the higher momentum 
transfers where no separation of CO and C2 is possible, we show in fig. |14| the structure 
functions A{q) and B[q). At the highest momentum transfers, the non-relativistic cal- 
culations with relativistic corrections may increasingly get doubtful, while the covariant 
calculations still could be valid. The comparison between theory and experiment, however, 
is not indicative of a real breakdown of the calculations which start from a non-relativistic 
frame work. The two-body terms that are poorly controlled — such as the p7r7 term — 
are mostly responsible for problems at large q. The calculations that use a soft pn'-j form 
factor do better than the ones using a hard one such as given by the VDM. 

While most of the calculations use the Hoehler nucleon form factors, the calculation of 



90| uses the dipole parameterization. The Hoehler Ggp, the most important nucleon form 
factor entering the calculation of the deuteron charge form factors, actually does quite 
well up to g ~ 9fm~^ in comparison to the precise Gep/Gmp ratios recently measured at 
JLAB by Jones et al. 0; the comparison between data and theory shown in this section 
thus should not be affected by shortcomings in the proton form factors employed. The 
exception might be the calculation of A{q) of ref. at the highest g's, where the use of 
the dipole parameterization does lead to too large an A{q). 

The success up to large q of the calculations that start from a basically non-relativistic 
frame work — supplemented by f ^/ relativistic corrections — may at first sight come as a 
surprise. Upon reflection, this may not be so astonishing, however, as the relation between 
q and the typical nucleon momenta probed is somewhat special for the A=2 system. For 
a heavy nucleus the typical momenta probed are of order q/2: when describing the form 
factor using a momentum space wave function, the biggest overlap between initial and 
final state occurs for \l/(— g/2)\l/(— g/2 -|- q). For the deuteron, the center-of-mass motion 
leads to the fact that this overlap is maximal at initial nucleon momenta ~ g/4; for the 
momentum transfers reached by experiment, g/4 does not yet correspond to genuinely 
relativistic momenta. 

A source of uncertainty in the prediction of the deuteron form factors resides in the 
contribution of AA components. Various groups have studied wave function components 
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involving the A or higher nucleon resonances ||164||-|T7T| using N — A transition potentials 



derived within OBE or boundary condition models. The calculation of the A components 
turn out to be difficult, and partly ambiguous, as the A^A^ A A transition is of very 
short range; the corresponding cut-off and vertex form factors are poorly known. As a 
consequence, there is little agreement between the calculations. While several approaches 
give AA probabilities in the 0.5-1% range, some give much larger probabilities. There 
is no agreement either on whether the contribution to e.g. A{q) is positive or negative 
at low q. It seems clear that the A A components have a larger effect on B{q) (near the 
diffraction minimum). It also emerges that the effects of A A components are significantly 
smaller than the ones of the standard two-body currents. Given the uncertainties of the 
latter, it is very difficult to learn from the comparison of theory and experiment something 
on the AA contribution. 

The calculations compared above to the experimental results obviously only represent 
a small fraction of the presently available theoretical predictions. Calculations of the 
deuteron form factors within different theoretical frame works are available: for more 
detailed consideration of alternative relativistic approaches such as the light front [p.72|| 



and equal-time approximation to the solution of Bethe-Salpeter equation ||173|| , we refer 
the reader to the original literature. 

A rather different framework has been employed by Buchmann, Yamauchi and Faessler 



1 74|| . These authors use the quark cluster model and the technique of the resonating- 
group method. The deuteron wave function is derived from a microscopic six-quark 
Hamiltonian which, in addition to the quark confining potential, includes a one-pion and 
one-gluon exchange potential between quarks. In this calculation the electromagnetic 
current operators are constructed directly on the quark level. 

The particular interest in this approach lies in the fact that, due to the application 
of the Pauli principle on the quark level, there are new quark interchange terms in both 
the one- and the two-body current matrix elements that are not present in a nucleon- 
level calculation. The occurrence of these terms can perhaps best be used to gauge the 



relevance of a six-quark (rather than two-nucleon) description of the deuteron. Fig. 15 
shows some of the results of this calculation. While it is probably not appropriate to 
quantitatively compare the one-body results to the data, the dashed curve shows the new 
exchange terms that appear in addition to the standard ones (dotted line). 

These additional terms clearly become important at the very large momentum trans- 
fers. Similar results have been obtained by Buchmann et al. for the Ml and C2 form 
factors, where the contribution of the quark exchange currents is smaller, however. 

The deuteron form factors in the past also have been discussed in terms of a very dif- 
ferent approach. The large momentum transfers achieved experimentally, and the smooth 
fall-off of A{q)., made it tempting to look at the asymptotic behaviour of the deuteron 
form factors in terms of the quark model ||175|| (the asymptotic behaviour in terms of 
nucleonic constituents [[L76| , |177|| curiously received much less attention). 

These predictions, starting from quark counting rules, were applied despite the fact 
that the condition for the applicability, q/6 ^ 3fm~^, where 3fm~^ is the typical mo- 
mentum of quarks and 6 is the number of (valence quark) constituents in the deuteron, is 
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Figure 15: 

quark cluster model predictions of Buchmann et al. 



Data for the deuteron monopole form factor compared to the 

Hi- 



Dotted curve: 



conventional pion exchange current, dashed curve: additional quark exchange 
current, solid: full calculation. 



far from true (12/6 ^ 3). These considerations on the asymptotic behaviour also contin- 
ued to be used despite the fact that the estimates for the 6-quark-contribution at the g's 
of interest were orders of magnitude below the data [|178|| . Today, with the experimental 
proof that the form factors display diffraction features (not g"" behaviour) and are neg- 
ative at the largest q (while the counting-rule predictions are positive as the normalizing 
constant is a probability), these quark counting guesstimates have fallen somewhat out 
of fashion. 



3.6 Electrodisintegration 

In addition to the elastic form factors, there is one transition form factor to a nearly 
discrete state of interest when discussing the form factors of few-body nuclei — the Ml 
transition to the singlet-S state of the deuteron right above disintegration threshold. This 
transition form factor has gained particular prominence as it has become one of the 
showcases for the role of two-body currents. We therefore include it in the discussion 
here. 

The first data for the Ml transition were measured by Rand et al. at 180° in conjunc- 
tion with a measurement of the deuteron magnetic form factor ||113|| . These measurements 
reached a momentum transfer g of 3.1 fm~^ and showed, despite the limited resolution, 
the enhancement due to the singlet-S state. These data, however, did not attract partic- 
ular attention at the time, as they appeared to agree with the then accepted lA-theory 
of Jankus. 



Only later on it was realized by Hockert et al. ||179|| that, similar to the A=3 magnetic 
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form factor, the transition form factor in lA should display a pronounced minimum near 
3.5/m~^ due to the destructive interference of the ^S-^S and ^D-^S transitions. This D-S 
transition had not been included in the Jankus calculation, as the deuteron ground state 
was assumed to be a pure S-state. The occurrence of this minimum leads to a pronounced 
disagreement of lA-calculations with the data, by a factor of 10. Only when allowing for 
two-body currents agreement with the data can be restored. As for A=3, the two-body 
currents give a basically opposite contribution to the S-D transition. The large effect 
of the two-body terms is essentially a consequence of the fact that electrodisintegration 
corresponds to an Ml isovector transition. 
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Figure 16: Inclusive cross section for d{e,e'), displaying the elastic peak 
{Ex = 0) and the enhancement at low E^; due to the transition to the singlet-S 
state p7| . 

Detailed measurements of deuteron electrodisintegration then were performed at low 
q at Mainz ||180| , |181|| who made an L/T-separation up to 2 fm~^. Measurements up to 



high q were done at Saclay by Bernheim et al. and Auffret et al. ||117| , |182|| who measured 



the cross section at a scattering angle of 155°. This experiment provided precise data up 
to momentum transfers of 5.2/m~^, thus covering fully the region where the lA-minimum 
and maximum of the cross section would occur. 

Fig. ITBI shows the enhancement of the cross section right at threshold due to the 



^S-state. In fig. the data are compared to an early calculation of Mathiot ||183|| . 



The lA calculation, accounting for the nucleonic contribution only, shows a pronounced 
interference minimum which does not occur in the data. The pionic two-body term gives 
the main contribution that brings theory closer to experiment. The p-exchange currents 
and isobar components produce non-negligible contributions as well, but the quantitative 
contribution was very sensitive to the cut-off masses assumed for the vertex form factors. 

When comparing data and calculations in general, some caution is advisable. While 
the data in most cases are averaged over the region of energy i?„p=0-3 MeV above 
disintegration threshold (0-10 MeV for the SLAG data), the calculations often consider 
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Figure 17: Cross section for deuteron electrodisintegration in the region of 
the lA minimum, averaged over the range E„p = 0-3 MeV [|117|| , compared to 
the calculation of Mathiot |183| (for E„p = 1.5 MeV). 



a fixed energy of 1.5 MeV above threshold. Given the cusp-shaped behaviour of the cross 
section near threshold (see fig. 16), this is not entirely equivalent and can, depending 
on momentum transfer, lead to appreciable differences. In particular, near the minimum 
at g = 5/m~^, the minimum for E„p = 1.5 MeV is much sharper than the one in the 
excitation-energy averaged data which includes contributions from multipolarities other 
than Ml which dominate at the larger Enp- 

Additional data have been provided by an experiment performed at Bates ||184| | and 
the NE4-experiment at SLAG ||185|| . This latter experiment, which reached the highest 
momentum transfers, was optimized to measure at 6 = 180° the elastic magnetic form 
factor, and had a fairly poor energy resolution (as much as 20 MeV). After unfolding 
of this resolution the cross section over the interval 0-10 MeV could be extracted. The 



complete set of data referred to above will be used in figs. |T8|-^ when comparing to theory. 

The calculation of the two-body contributions has made considerable progress since 
the early calculations cited above. In particular, it is now possible to derive the vr and 
the p-like exchange currents consistently with the nucleon-nucleon potential employed for 
the calculation of the nucleonic part. Also, the relativistic effects of order v^/c^ can be 
fully included, thus removing the earlier ambiguity of whether to use Ge or Fi as vertex 
form factors. The isobar contribution continues to be somewhat uncertain, as it can not 
be fixed via the N-N potential and current conservation. 



Fig. shows, for a modern calculation ||186|| , the different contributions to the cross 



section, averaged over the region Enp= 0-3 MeV. Glearly, the two-body terms lead to 
the biggest change relative to lA in the region 2-6 fm~^. Isobar currents (IG) somewhat 
increase the cross section without changing much the overall g-dependence. Relativistic 
effects (RG) become very important at the largest momentum transfers. 
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Figure 18: Prediction for d(e,e') cross sections up to large momentum transfer 
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Figure 19: World data for d(e,e') integrated over E„p-regions as indicated. 
The cross sections are compared to calculations including MEC of Mathiot 
181 and the group of H. Arenhovel ||186||, for two different N-N potentials. 



As shown by fig. |T9| the calculated cross sections are reasonably sensitive to the N-N 
potential employed. This is mainly due to the different predictions for the D-state, where 
the energy-dependent Bonn potential produces quite different results as compared to e.g. 
the Paris potential. One also should note that the calculated cross sections are not too 
sensitive to the nucleon form factors employed ||150|| as disintegration at threshold is a 
magnetic transition, so Gen is not important. 

The large momentum transfers reached, together with the ambiguities of the calcula- 
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tions of the two-body terms occurring when deahng with inter-nucleon distances that are 
comparable to the nucleon size, have led to alternative attempts to calculate the disinte- 
gration cross section. The results of some of these approaches are displayed in fig. 120. Lu 



and Cheng [|187|| use a hybrid model ||188|| where the deuteron is divided into two distinct 
regions: an exterior one described in terms of baryons, and an interior one described in 
terms of six-quark configurations. The matching radius chosen was near Ifm. This model 
predicts a second maximum of the cross section at g ~ 7fm~^, a maximum taken as a 
signature for a six-quark contribution. This maximum, however, was not seen in the data. 
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Figure 20: Cross sections for d(e,e') predicted using quark-based models 
189| , |187|| , or a chiral Lagrangian ||190 |. 



Quark-based calculations have been published by Yamauchi et al. ||189|| and Chemtob et 
al. ||191|| . In the calculation of Yamauchi et al. the two-nucleon system is described using 
the resonating-group method. The quark-quark interaction includes quark interchange 
effects, with gluon-exchange corrections, and pion-exchange effects. While the calculation 
is quite successful at large momentum transfer, it predicts too pronounced a minimum at 
q = 5fm~^, which would only partly be filled in when averaging over the experimental 
E„p-range. 

Smejkal et al. [ |19CI|| construct a Lagrangian in terms of N, vr, p and Ai. The exchange 
terms include vr and p exchange and Ai — n exchange. The resulting cross sections are 



also shown in fig. |20. 

The comparison to the data shows that the calculations in terms of nucleons with 
added two-body terms consistently derived from the N-N interaction are quite successful 
in explaining the data. The calculations in terms of quark constituents for the time being 
do not correctly reproduce the data. Part of the success of the conventional approach 
may again be due to the special nature of the deuteron mentioned already above. The 
typical internal momenta probed in the deuteron wave function are of order g/4, i.e. less 
than 2 fm~^ for electrodisintegration; at these momenta, the conventional description in 
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terms of nucleons and mesons still can be expected to work well, while quark-model based 
approaches could be expected to be valid at much higher g's only. 



4 3-body nuclei 
4.1 Introduction 

Three-nucleon elastic form factors are particularly useful for assessing our understanding 
of the few-nucleon systems and our abihty to quantitatively predict nuclear properties 
based on the knowledge of the nucleon-nucleon interaction. The wave functions and 
form factors can reliably be calculated starting from modern nucleon-nucleon interactions 
known from the two-nucleon system; two-body exchange currents can be taken into ac- 
count in a way that is largely consistent with the nucleon-nucleon interaction. Apart from 
a much higher nucleon density, the main difference between the two- and three-nucleon 
systems is related to the presence of the third nucleon, which can lead to the appearance of 
three-nucleon forces. The most sensitive observables to test these theoretical predictions 
are the electromagnetic elastic form factors. The beak-up cross section at low energy loss, 
which in analogy with the deuteron potentially could also provide interesting information, 
has hardly been exploited up to now ||192| . 



For momentum transfers below 5fm the one-body contribution (impulse approx- 
imation) is now considered to be quantitatively under control. It is a well known fact, 
however, that the impulse approximation alone can not explain the magnetic or the charge 
form factors. Two-body currents as well as relativistic effects have to be added in order 
to get satisfactory agreement with the data. 

All trinucleon form factors are accurately known up to momentum transfers q of the 
order of 5.5 fm~^, transfers large enough to cover the region of the secondary diffraction 
maxima. These data are now adequate for an extensive comparison with the present 
theoretical calculations performed within the framework of nucleonic constituents and 
two-body currents. Up to the momentum transfers where the ^He and form factors 
reach similar values of q, this comparison can be extended to the combinations of isospin 
T=0 and T=l, thus allowing a direct confrontation with form factors for the 2- and 4- 
body systems. These include the deuteron charge and magnetic form factors which are 
both pure AT = transitions, and deuteron electrodisintegration at threshold, which is 
a AT = 1 transition. 

The ^He and magnetic form factors are of particular interest due to the destructive 
interference that occurs at large q, where the S-S and the S-D terms partly cancel pB|. 



This leads to a pronounced shift in q of the predicted interference minima. The A=3 
magnetic form factors are also of special interest when studying the effects of two-body 
currents. Already at g = their contribution is sizeable: the magnetic moments in lA for 
^Hi^He) are 2.57 (-1.76) magnetons while the moments including two-body effects 
amount to 2.98 (-2.09). 
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4.2 Electron scattering data 

The and ^He nuclei present particular challenges for the experimentalist. For both 
nuclei targets of adequate thickness are difficult to produce and use, particularly of course 
for the radioactive ^H, which we discuss first 

The pioneering experiment on was performed by Collard et al. at Stanford ||193|| . 



For both and ^He high-pressure (lOObar) sealed gas targets in the form of cylinders 
of 2cm diameter and 20cm length were used. The target contained as much as 25 '000 
Curies. The and ^He data could be measured quasi- simultaneously, thus ensuring 
low systematic errors in the relative cross sections. Data were taken up to a momentum 
transfer of 2.8/m~^ for both the charge and the magnetic form factors, so the data did 
not yet reach the diffraction feature located at larger q. 

An experiment using a somewhat different technology for the gas-target was carried 



out more than 20 years later by Beck et al. ||194|| using the Bates accelerator. Vertical 
cylinders of ~10cm diameter were used as target. The was stored in a Uranium oven, 
and transferred to the target which was kept at 45K temperature and a pressure of 15 bar, 
yielding a target thickness of ~ bOmg/cm'^. This target contained a rather large quantity 
of ^H, 140 kC. Contrary to the target of ref. ||193|| the target-system did require handling 



of the radioactive gas in the experimental hall. This experiment, which also reached 
momentum transfers of ~ 2.8 fm~^, did achieve higher accuracy as the luminosity was 
considerably higher than the one utilized in ref. ||193|| . The use of not too far from the 
point of liquification did, however, introduce some uncertainty on the target thickness as 
the thermodynamical properties of have to be extrapolated from Hydrogen. Again, 
and ^He data were measured in the same experiment. This experiment also measured 



quasi-elastic scattering data ||195 |. 



For the Saclay experiment ||196| , |197|| a very different technology was utilized. Amroun 



et al. used a target containing liquid tritium, as only this option provided the luminosity 
required to reach momentum transfers beyond the expected diffraction features in Fch 
and Fjn. For reasons of safety, a sealed cryogenic system was employed that allowed to 
liquefy without any radioactive-gas handling, except for the initial fill and final gas 
recuperation carried out in a specialized laboratory. Again for safety reasons, the tritium 
target was enclosed in three further sealed containers with thin windows for the electrons, 
and a massive container closed during transportation of the target. The general layout of 
the target is shown in figure PT| . 

The ^if-target proper consisted of a stainless steel cylinder, 40mm long, 10mm in 
diameter with hemispherical end caps located outside the acceptance of the spectrometer. 
The target was in permanent connection with an expansion- vessel of ISOcm^ volume. This 
vessel was equipped with a heater, in order to keep it at T ~ 300K while the target was 
cold. When the entire system was at room temperature, the tritium pressure was 23 bars. 

This combination of cold target and warm storage vessel had a number of outstanding 
benefits, such as high target density (liquid ^H), low operating pressure while the target 
was placed in the high- intensity electron beam (3 bars), very high efficiency in the use 
of the (98% in the cold target cell), well known density of liquid ^H, and thin target 
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Figure 21: Saclay tritium target. The volume between containers C3 and C4 
is filled with gaseous ^He for a continuous check of integrity of C3 and C4. 
The volume C4 could retain the tritium at pressures below 1 bar. 

windows. Only lOkC of tritium were needed to reach the high target thickness. 

Use of a liquid target with minimized amounts of involves a large amount of heavy 
nuclei (Cu, Fe) near the beam. A uranium collimator (diameter 8 mm) upstream of the 
target and outside the spectrometer acceptance eliminated potential stray electrons. A 
secondary emission monitor upstream of the target verified that no beam-halo was present. 
A split secondary emission monitor ensured stabilization of the beam position to 0.1 mm 
accuracy. 

For the extreme forward- and backward-angles (25° and 155") the end caps of the 
target were not entirely outside the acceptance of the spectrometer. Insisting on this 
condition would have led to an undue increase of target length {i.e. quantity of ^H). For 
the extreme angle of 155" — which is of great interest to the separation of Fch and Fm 
— a special collimator placed close to the target ehminated all contributions from the 
end caps. 

At the nominal temperature of 21. 7K the target density was 0.2710±0.0014 gjcm?. 
This led to a thickness of ~ 1 gjcrr? at 6* = 155°. With this large target thickness, 
data could be taken at momentum transfers where the cross section is very small. As a 
consequence, the region of the expected diffraction minimum and secondary maximum 
could be covered, and the experiment reached a maximum q of 5.7 fm~^. 
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For ^He the life of the experimentahst is somewhat easier, as no precautions due to 
radioactive gas are required; achieving an adequate target thickness still is a challenge. 

Besides the experiments mentioned already above |193| , |194| several groups have con- 
tributed to the ^He data base. An experiment performed at Mainz by Ottermann et al. 
198|| provided data up to g = 1.9/m~^, lower-g data were measured by Szalata et al. [[199 



and Dunn et al. 200 



The region of the diffraction minimum and secondary maximum of the charge form 
factor was first reached by the experiment of McCarthy et al. ||201| , |202|| who, contrary 
to all other experiments, used a liquid target cooled to IK using superfluid Helium. The 
diffraction feature in the magnetic form factor was reached in the experiment of Cavedon 
et al. ra 



al. 



The data on the charge form factor going to the largest q were taken by Arnold et 
[P04|| in a coincidence experiment performed at SLAC. Due to the cryogenic target 
(42cm of ^He at 20K and 50 bar) a large target thickness was achieved, ~ ^g/cm?. In 
combination with the high energy of £'e=15 GeV used — the cross section at a given q is 
roughly proportional to — very small form factors (F ~ 3 • 10~^) could be measured. 
This allowed Arnold et al. to reach a maximum momentum transfer of 10/m~^, and cover 
the region of the second diffraction zero. 

The form factors presently available, as extracted directly from the various experi- 
ments, are shown in fig. ^ 
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Figure 22: Charge and magnetic form factors of ^He. 



An experiment performed at Bates ||205|| has produced two points at higher q for 
magnetic scattering, the data have, however, not yet been published. 
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4.3 Experimental form factors 

The world set of cross sections measured for ^He and have been analyzed by Amroun 



et al. ||197|| in the same way as described in sect. p.3| for the deuteron. The various form 
factors were parameterized using the SOG parameterization. From these form factors, 
the cross sections were calculated and the parameters fitted to the data. The separation 
of longitudinal and transverse form factors thereby was automatically achieved. 

In order to analyze the electron scattering data it was assumed that the cross section 
can be expressed in the plane wave Born approximation (PWBA) as: 



CMott- 



Fl{q) + 



w 

2M2 \ 2q 



where the Mott cross section contains the usual recoil and factors, /i is the ^H{^He) 
magnetic moment, rj is given by (1 + g^/4M|.). M (Mj-) is the nucleon (three-nucleon) 
mass and q (q) is the four (three)-momentum transfer. 

As the incident electron wave is distorted by the Coulomb field of the nucleus, the 
PWBA is not entirely valid. The DWBA form factors F^f^ j^^y{E 
related to the Born form factors F'^^^q) and F^{q) by: 



and Fl^jj^{E, 



are 



p2 

^ch 



Fm,DwiF,(^)fm{E, 



where fc and fm are Coulomb correction factors. Unfolding of these corrections allows to 
restore the PWBA formula and to infer the PWBA form factors that can be compared 
to theoretical predictions and will be used below. 
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Figure 23: Best fit charge form factors (curves with error band) and form 
factors derived using standard Rosenbluth plots in individual regions of q. 



In order to unfold the Coulomb corrections from the experimental cross sections ap- 
proximate charge and magnetization densities were used to compute the cross sections as 
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a function of the actual energies and angles. This was done using a phase shift code for 



charge scattering ||206|| , and the DWBA code HADES [P07|| for magnetic scattering. The 
Coulomb correction factors fc, fm, obtained from the ratio of the exact to the PWBA 
cross sections, were found to be insensitive to the exact densities employed. 

Amroun et al. have also performed standard Rosenbluth separations of the data. As 
forward- and backward-angle data are not exactly at the same momentum transfer (al- 
though close as by now data at many angles are available), they have used the fit to 
shift the data to the g-value desired. At the g-value of interest the ratio aexplo'fit of the 
neighbouring experimental g-value was used. The resulting Fch and Fm agree with the 
ones from the global fit. The error bars are less precise, however, as for a determination at 
a given q the global fit de facto uses the values of all data in the neighbouring g-interval; 
the width of this interval is given by ~ 1/Rmax — 0.2/m~^. 

In the paper of Amroun et al. ||197|| the best fit parameters for and '^He are given 
together with the uncertainties of the form factors, so we dont need to repeat them here. 
(Note that eq.(l) contains a misprint, the numerical factor in the exponent should be 1/4 
and not 1/2). 

In fig. ^ the comparison between the global fit (curves with error band) and the form 
factors determined from Rosenbluth separations in individual intervals of q is shown. 
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On occasion, it is useful to know the rms-radii of the A=3 nuclei. They are listed in 



table |4.3| , where the error bars given include both statistical and systematic uncertainties. 
For the magnetic form factors, the "radii" given do not really correspond to a genuine 
radius; rather, they simply reflect the coefficient that occurs in the g = expansion of 
the magnetic form factor F(g) = 1 — g^r^^g/6 + .. 



4.4 Comparison to theory 



As mentioned in sect. |2.3| a number of techniques are available to calculate the wave 
function starting from the (non-relativistic) Schrodinger equation and the known N-N 
interactions: Faddeev calculations in both coordinate- and momentum space B^, IBB, p7|. 



20^ , |209| , |210| , |2 1 1| , pi2|| , the Correlated Hyperspherical Harmonics approach [|7T], |213|| and 
Greens-function Monte Carlo methods pi4 . Relativistic calculations are also becoming 
available ||215| , pi6| , |49| , |217|| . Calculations of the A=3 nuclei have been performed using the 
Blankenbecler-Sugar formalism |^ and the relativistic spectator (Gross) equation ||217|| . 
These relativistic calculations are of rather recent date, however, and the corresponding 
electromagnetic form factors have not yet been published. 
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In the following we first compare to the results of the recent calculation of Marcucci et 



al. [21S]. These authors use the trinucleon wave functions calculated by Kievsky et al. |71 



with the Pair-correlated Hyperspherical Harmonics (PHH) method. In this PHH calcula- 
tion the AV18 N-N interaction, which includes charge-symmetry and charge-independence 
breaking terms and gives a very good fit to N-N scattering data, is used. 

The AV18 force is complemented by the UIX three-body force which has a long-range 
part given by 27r-exchange. At short range a phenomenological repulsive term is added 
to simulate dispersive effects. This part is fitted to the binding energies of the A=3, 00 
systems. 




q (fm 



Figure 24: Contributions of different two-body charge terms, split into T=0 
and T=l components pi8||. The one-body part is indicated by 1-N. 



The two-body operators are calculated as discussed in sect. p.4| , i.e. they are con- 
sistently derived, as far as presently possible, from the N-N interaction employed. For 
the two-body currents the vr- and p-like terms (PS and V in fig. ^4|) come from AVIS. 
The mo del- dependent terms, 7rp7 and Ci;p7, are calculated using known coupling constants 
and VDM form factors. The cut-off parameters needed at short range have been chosen 
to give a good description of the deuteron magnetic form factor B{q) at high q. The 
three-body currents resulting from the 27r-exchange in the 3BF employed is accounted for 
as well. For the two-body charge operators, which do depend on how the non-nucleonic 
degrees of freedom {NN, N*) have been eliminated from the wave function, the authors 
follow as much as possible the derivation of the corresponding vr- and p-like terms in the 
two-body currents. Again the npj and urpj terms have to be added separately. In fig. ^ 
the different contributions of the two-body charge terms are shown, split into their T=0 
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and T=l components. It is obvious from the comparison to the one-body curve that 
two-body contributions will have an important effect upon the form factor in the region 
of the diffraction minimum and maximum. 




Figure 25: 
et al. [BTS 



A=3 charge form factors together with calculation of Marcucci 
The lA result (dashed) , and the full result including two-body 



terms, is given. 



Marcucci et al. have also explored the role of A-components in the A=3 ground state. 
Often these degrees of freedom are included via the two-body currents, but a different 
picture might result once they are explicitly allowed for in the wave function, such as 
done previously by the group of Sauer ||219| , |220|| . Marcucci et al. do not carry out a full 
coupled channel N-A calculation. Rather, they use the Transition Correlator Operator 
(TOO) method ||221|| , which has proven to be useful in earlier variational calculations. In 
this case, they use the AV28Q N-N interaction which allows for both nucleons and Deltas. 



51 




Figure 26: A=3 magnetic form factors together with calculation of Marcucci 



et al. ||218|| . The lA result, the full result including two-body currents, and the 
result obtained when explicitly including A's are given. 



In figs. |25|j26| we show the comparison between the experimental results of the A=3 
form factors and the predictions of Marcucci et al.. In each case we show both the lA- 
prediction, and the full calculation including two-body and three-body currents. For the 
case of the magnetic form factors, the results obtained by explicitly including the A in 
the wave function are also shown. 

One observes that for all form factors the effect of two-body currents is quite important. 
In general, the full calculation does a very good job in explaining the data. The effect of 
explicitly including the A is rather small, significantly smaller than what was found by 
Strueve et al. ||220|| , partly because different coupling constants were employed. 

In figs. we compare the experimental data with further calculations of refs.| 219 , 

220| , P22| , |223| , |224|| , which also are based on nucleonic degrees of freedom. 

Sauer and collaborators ||219| , f^20| have solved the coupled-channel N-A Faddeev equa- 
tions in momentum space, and used the Paris NN potential modified ||219|| to include A- 
isobar excitations via vr- and p-meson exchanges. The presence of the A thus accounts for 
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the most important part of the three-body force. In this calculation tt-, p-, and p77r-meson 
exchange current contributions are included as well as the one-body Darwin-Foldy term 
and spin-orbit relativistic corrections. As for the calculation of ref. ||218|| boost corrections 
are not yet included. 



Schiavilla and collaborators |P22| , 223] based their (older) calculation on a variational 



three-body wave function computed using the Argonne V14 interaction and the Urbana- 
VII three-nucleon force. In this calculation the vr-like and p-like meson exchange cur- 
rents were derived consistently from the NN interaction used. The authors also included 
the one-body relativistic corrections mentioned above as well as a;-exchange and a;7p- 
contributions. 

Golak et al. ||224|| in their recent calculation solve the Faddeev equations in momentum 
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space, using the AV18 N-N potential. The two-body currents, which up to now are only 
calculated for the magnetic form factors, are treated according to ref. 
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For all calculations the two-body terms significantly displace the diffraction minimum 
of the charge form factors to lower q values and increase the height of the secondary max- 
imum by more than a factor of two. The calculations achieve a similarly good agreement 
with ^H. This is not the case for '^He where ref. | |220| | slightly underestimates the charge 
form factor, whereas the calculation of ref. ||222|| almost perfectly agrees with it. The 
difference could be due to the model dependence of the one-body term, or due to different 
two-body contributions. Indeed, for the form factor, the calculations of refs. ||220|| and 
223|| predict the one-body minimum at 3.9 and 4.15/m~^ respectively, while the predic- 



tions for the full calculation perfectly agree. Additional differences might come from the 
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different parameterizations of tlie nucleoli form factors used, particularly for the poorly 
known Gen- We do not understand why the form factors are better predicted than the 
ones for ^He. 

For the magnetic form factors, where the effect of the three-body force is also quite 
small, the effect of the two-body currents is larger than for the charge form factor, and they 
improve considerably the agreement with experiment. However, none of the calculations 
perfectly accounts for the shape of the '^He form factor: the prediction of ref. | 219 | explains 
its shape below the diffraction minimum, but fails above, whereas ref. [p22|| shows the 
opposite behaviour. The latter calculation agrees with the form factor, while this is 



not the case for ref. |pi9 |. 



There are several differences between these two calculations. As for the charge form 
factors they differently predict the one-body contribution. The lA predictions for the 
location of the diffraction minima of are 3.5 and 3.7 fm^^ for the calculations of 
Sauer and Schiavilla, respectively. The differences in the exchange current contribution 
are even larger. 

While qualitatively the role of two-body currents is understood, the situation for the 
magnetic form factors in particular is not yet entirely satisfactory when it comes to a 
quantitative comparison. Even at low q, some differences between experiment and theory 
remain. While experimentally the i^He) magnetic moments amount to 2.98 (-2.13) 
/iAT, the theoretical predictions are e.g. 2.98 (-2.09) for Vis and 2.87(-2.06) for the Paris 
interaction. 

Fig. E9| shows, for completeness, the isospin-separated A = 3 charge form factors 



together with the theoretical predictions of Strueve et al. and Marcucci et al. ||22CI| , |218 
The main observation from these comparisons is that the T=0 charge form factor is 
remarkably well predicted by both calculations. This agreement is quite unexpected since 
it is reached after adding "model-dependent" two-body contributions, which lead to a 
sizeable change. For the T=l charge form factor we observe larger differences between 
the theoretical predictions. Moreover, none of them perfectly describes the data. Since 
the agreement for the T=0 component has been found to be excellent for all of them, we 
conclude that improvement in the description of the trinucleon requires improvement of 
the calculation of its isovector charge component. 

For the T=0,1 magnetic form factors the calculations are closer to the data, so we do 
not separately show them. 



4.5 A=2,3 comparison 

In order to understand better whether eventual shortcomings in our understanding of 
the 3N-systems are only related to the two-body input, or whether our understanding of 
genuine three-body properties is lacking, it is of interest to directly compare the 2- and 
3-body form factors. 

To achieve this, one can extract from the ^He and ^if form factors the corresponding 
isoscalar- and isovector quantities, where the former then is directly comparable to the 
deuteron form factor. This comparison becomes particularly interesting now that the 
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Figure 29: T=0 and T=l charge form factors of the A=3 systems, compared 
to the calculations of Marcucci et al. |pl8|| (solid) and Strueve et al. ||220| 
(dashed). 



diffraction structure of the A=2 CO charge form factor is well determined. 

The isoscalar Ft=q and isovector Ft=i components of the A=3 charge form factors are 
defined by the following expressions: 

Fo(g) = \ [2FHe{q) + Fuiq)] F,{q) = 1 [2F^,(g) - F^(g)] 

and similarly for the magnetic form factors: 

1 1 

Foiq) = 2 [f^HeFneiq) + iJ^HFniq)] i^i(g) = - [iJ^HeFneiq) - f^HFHiq)] 

where /ih sjid fine are the A=3 magnetic moments and all ^He and form factors 
are normalized to 1 at g = 0. The T=0,1 form factors can easily be calculated from the 



parameters of the He and H form factors |197|. The effect of the isospin impurity, due to 



the Coulomb interaction in ^He, has been neglected. Indeed, its effect on the ^He charge 
form factor is very small, < 1% in the region of the secondary diffraction maximum and 
even smaller for lower q values. 



Henning et al. ||155|| have made a detailed comparison between some of the global 
properties of the A=2,3 charge form factors such as the position of the diffraction minima, 
the position of the maxima and the amplitude of the form factors in the maxima. They 
have done so by calculating the A=2,3 form factors for a number of N-N interactions, 
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Figure 30: Correlation between the position of the diffraction minima in the 
A=2 and A=3 isoscalar charge form factors (top) and the amphtudes of the 
form factors in the diffraction maxima (bottom). The experimental point is 
indicated by the shaded area (exp). 



with and without inclusion of the two-body currents. The potentials employed were the 
venerable RSC potential, the Paris potential and various OBE Bonn potentials OBEPQ, 
OBEPA,OBEPB and OBEPC. 

It turns out that the above-mentioned observables display a remarkable correlation 
between the A=2 and A=3 systems. There is a basically linear relationship between the 
A=2 and A=3 quantities. This is demonstrated in fig. ^ which gives the amplitude 
of the A=3 isoscalar charge form factor in the diffraction maximum as a function of 
the same quantity for the deuteron. Not only do the values for the different potentials 
(to which we have added the AV18 potential) define a basically linear correlation, the 
points with/without two-body currents also lie on the same line. Within the line, there 
is a correlation between the strength of the tensor interaction and the position of the 
minimum: N-N potentials yielding a larger D-state produce a minimum at lower q. The 
lower part of the figure gives the corresponding relation between the positions of the 
diffraction zeroes. 

The inclusion of the AV18 points shows that this line is rather general; it is not simply 
a consequence of the particular way the two-body contributions are calculated. There is 



some qualitative understanding for the occurrence of this correlation ||155|| : if the A=3 



57 



nuclei could be described in a quasi-deuteron model, such a correlation would be expected. 
A quasi-deuteron model is obviously much too primitive for the A=3 system, but the three- 
nucleon wave functions obtained from the Faddeev equations do approximately show the 
symmetry upon interchange of the Jacobi momenta p and \/3/2q which would result from 
the spectator model. 

The non-understood feature of fig. ^ lies in the fact that the line describing all the- 
oretical calculations does not pass through the experimental point indicated by "exp", 
particularly for the amplitudes in the diffraction maxima. From the comparison with ex- 
periment it would seem that the existing range of non-relativistic N-N potentials together 
with the present understanding of two-body currents is not able to account simultaneously 
for the properties of the A = 2,3 systems. 

Besides the difficulties with the ^He magnetic form factor, the discrepancy shown in 



fig. 30 is one of the major open problems, and potentially points to the need to go beyond 



the "nuclear standard model" 



5 Helium-4 

The four-body system is interesting as the a-particle is tightly bound, and the nucleons 
are dominantly in a spatially symmetric S-state configuration. As a consequence, a high 
central density (the highest in nuclei) is found in ^He. 

The data base for the '^He charge form factor is reasonably complete, and the form 
factor measurements extend to unusually large q. The pioneering experiment was per- 



formed by R. Frosch et al. at Stanford | 226 | using a liquid helium target. This experiment 
already reached the diffraction maximum at 4/m~^. 

More precise data at medium momentum transfers were provided by an experiment 
carried out also at the Stanford HEPL by McCarthy et al. ||202|| ; this experiment focused 
on the ^He form factors, the '^He was measured mainly as a check. The separation of 
longitudinal and transverse contributions, necessary for ^He, was performed for ^ife as 
well, in order to certify that indeed a transverse form factor compatible with zero was 
found for '^He. As the systematic errors are very different for the forward angle (high 
energy) and backward angle (low energy) data, this provided a very stringent and useful 
test of the quality of the data. 

Data at very low q are also available from the experiment of Erich et al. [p271| performed 
at the Darmstadt electron accelerator. This experiment used a mixed H/He gas target 
in order to obtain a continuous normalization of the helium cross sections relative to the 
proton. Low-g data were also measured by Ottermann et al. [|198|] at Mainz, using a 14bar 
gas target, coUimated to exclude the contribution of the target windows. 

The highest momentum transfers were achieved by an experiment carried out at SLAC 



by Arnold et al. \\204\\ . In this experiment, which used electron energies up to 15GeV, 
a long gas target at 20K and 50 bar was used to obtain a high target thickness and 
large luminosity. Due to the high electron energy and the correspondingly poor energy 
resolution, the scattered electron and recoil helium had to be detected in coincidence 
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using the 20GeV/c and 8GeV/c-spectrometers. The measurements reached a momentum 
transfer of where the count rates dropped to 1 event/week; the steepness of the 

fall-off of the form factor at the highest q is indicative of the presence of a second diffraction 
zero. 

For ^ife, as for the other few-body nuclei discussed in this review, we have used the 
world data to determine the experimental form factor. The treatment of the data and 
their uncertainties is similar to what was discussed for the case of the A=2,3 systems. 
The resulting form factors are shown below in the figures. 

Several groups have calculated the ground-state wave function of the 4-body system. 
As for the A=3 nuclei, '^He is underbound by several MeV in calculations that only 
include two-body forces. The contribution of a three-body force must be added. As a 
matter of fact, the ^He binding energy is often used to fix the scale factor of this 3-body 
force. 

Glockle and collaborators [ [228| , |229 | have employed the Faddeev-Yakubovsky technique 
using a selection of N-N interactions, Viviani et al. ||213|| have employed the Correlated 
Hyperspherical Harmonics (CHH) approach together with the AV14 and AV14+TNI8 
interactions. As far as ground-state properties such as binding energies go, these calcula- 
tions compare well to the Variational Monte Carlo (VMC) results obtained by Arriaga et 



al. [230] and the Greens- Function Monte Carlo results of Carlson and Schiavilla | 231 , 232 |. 



Further results for the ^ife charge form factor are available from Schiavilla et al. [ p23|| 
based on a VMC calculation using the Argonne V14 potential supplemented by the Urbana 
VII three-nucleon interaction. The most recent result has been obtained by Schiavilla and 
Carlson using the VIS potential supplemented by the Urbana IX three-nucleon interaction 
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Figure 31: Experimental form factor of ^i?e together with lA (dashed) and 
full calculation of Carlson and Schiavilla ||233|| performed using the Argonne 
V18 N-N potential together with the Urbana IX three-body force. 
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In fig. ^ we show the results of the latest calculation together with the lA-prediction. 
The figure shows that modern calculations do extremely well in reproducing the experi- 
mental data. The effect of the two-body terms is very important in achieving this good 
agreement, particularly for the height of the diffraction maximum at 4/m~^. This is the 
case despite the less than satisfactory state of the theory of isoscalar two-body contribu- 
tions for the charge form factor. When comparing the isoscalar charge form factors of 
A=3 [|197|| and A=4 one notes that in both cases theory somewhat overpredicts the form 
factor in the diffraction maximum. 

The breakdown of the different contributions to the two-body terms is given in fig. 12 of 
ref . p23|| . The individual contributions are quite similar to what we have shown in fig. E| 



for the isoscalar A=3 form factor. The vr-like (PS) and p-like terms (V) give the largest 
contribution at high momentum transfer. The single-nucleon relativistic correction, the 
Darwin-Foldy (DF) and Spin-Orbit (SO) terms are also significant at large g, and help to 
get large enough a (absolute value of the) form factor in the diffraction maximum. The 
contribution of the most model-dependent two-body terms fortunately is relatively small. 
In order to show the progress of the understanding of the '^He form factor, fig. ^ 



displays the result of an early calculation by Gari et al. [ p34|| . This calculation was 
performed using the exp(S)-approach, including up to 4-particle correlations. The curve 
shown was calculated using the RSC N-N interaction together with a phenomenological 
three-body force. The two-body contributions were calculated according to the approach 
by Gari and Hyuga ||148|| . The main discrepancy to the data occurred in the region of the 



diffraction maximum. The comparison to modern calculations (fig. |3T|) shows that both 
the lA contribution and the two-body contributions have increased by nearly a factor 
of 2, due to both a better ground-state wave function and a more consistent calculation 
of the two-body processes. In fig. |3^ we also show the results of an older calculation of 
Schiavilla et al. ||223|| obtained using the Argonne V14 plus Urbana-VII TNI interactions. 
A similar observation can be made for this calculation. 



The data have also been used by Sick |p35|| to determine an accurate charge radius, 
with the goal to compare it to a precise measurement from the 2p-2s transition in muonic 
helium by Carboni et al. | |236| |. For ^He such an analysis can be done by including 
information from other observables on the large-radius behaviour of the proton radial 
wave function. This large-r behaviour has a considerable influence on the rms-radius, but 
is not very well determined by electron scattering as at small g, where its effect shows up, 
the experimentally determined signal, the difference of the form factor to one, becomes 
small and affected by systematic errors. 

For '^He the shape of the proton radial wave function is calculable knowing the proton 
separation energy. The absolute normalization is known from a dispersion-relation anal- 
ysis [p371| of the world-data on proton-helium elastic scattering. This constraint helps to 
extract from the (e,e) data a very precise charge rms-radius, 1.676 ± 0.008 /m, which 
agrees with the value measured using muonic atoms. 



The calculation corresponding to fig. ^of Schiavilla et al. | [23ij | gives a charge radius 
of ^ife of 1.609 fm. Once the contribution of two-body currents is added, this increases 
to 1.634 fm. In comparison with the experimental result, this is a bit low. 
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Figure 32: Experimental '^He form factor together with early results from the 
exp(S)-calculation of Gari et al. ||234|| (solid curve) and the VMC calculation 
of Schiavilla et al. ||223|| performed using the V14 interaction (dashed). 



Overall, the understanding of the ^ife form factors is quite satisfactory. This shows 
that the calculations of the isoscalar two-body contributions are quite successful even in 
dense systems, although the theoretical understanding of these terms still leaves to be 
desired. Other open questions on the structure of ^He, such as the size of the D-state 
contribution or the role played by vs "^Si pairing, are better investigated by studying 
breakup observables. 



6 A>4 



The region of "light nuclei" up to a few years ago certainly did not go beyond mass 
number A=4: for heavier nuclei calculations of accuracy comparable to A<4 could not 
be performed. In this situation, one had to take recourse to more approximate, and much 
less fundamental, approaches such as the shell model, cluster models etc. Some of these 
calculations have reached quite a high standard (see e.g'. ||238|| ), but these approaches do 
not provide the direct link to the underlying N-N interaction. 

This situation has changed during the last years. Today, Greens- Function Monte Carlo 



(GFMC) calculations can be performed for nuclei up to mass A=9 |]64| , [76| , and Variational 
Monte Carlo (VMC) calculations have been performed up to A=16 ^ |25D[. A 
recent exp(S)-calculation also reaches A=16 [p41|| . In these calculations, the nuclear wave 
function can be obtained starting from the N-N interaction, with no additional model 
assumptions. Despite the daunting numerical effort required and some approximations 
that do have to be made, the results seem to indicate that a fairly accurate understanding 
is being reached for these not-so-light nuclei. 

It is of great interest to study the p-shell nuclei in terms of "exact" methods. For the 
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heavier p-shell nuclei, at least, the shell model already provides a reasonable description; 
in this case the comparison between the "exact" calculations and the model-calculations 
allows us to assess the accuracy and potential failures of the latter, and herewith lets us 
gain most valuable insight needed when dealing with even heavier nuclei where "exact" 
calculations are out of question for quite some time to come. 

GFMC calculations for nuclei up to A=9 have been performed by Pudliner et al. 
These calculations employed the AV18 N-N interaction together with the Urbana 



64, 76 



IX three-body force. For the binding energies, these calculations reach accuracies of the 
order of IMeV out of AOMeV. The calculations show that the N-N correlations continue 
to play an important role even for the heavier nuclei; the two-body densities for these 
heavier nuclei actually at N-N distances < 2 fm are quite similar to the ones found for 
the deuteron. 
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Figure 33: ^Li form factors calculated by VMC by Wiringa and Schiavilla (o) 
242| together with data for C0+C2 l24|, HI and Ml i24|, HI, (+). 



Only for some of the p-shell nuclei and selected states have the electromagnetic form 
factors been calculated. In these calculations the two-body currents have been included 

as far as possible consistently with the underlying 
In figs. p3],p4| we show two examples. For the ^Li charge form factor 



in the way discussed in sect. 
N-N interaction. 

the results calculated using VMC agree quite well with the data of G.C. Li et al. [p43 
measured at Stanford, and Bumiller et al. [M^] and Suelzle et al. | 248 | at the lower 
momentum transfers. The calculations include the contributions from two-body currents 
which have an effect only near the diffraction minimum where they shift it to somewhat 
smaller q. Both the data and the calculation correspond to the sum of the CO and 
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C2 multipolarities {^Li has spin one) which have not been separated. The agreement 
between calculated (2.59±0.02/m) and experimental (2.56±0.05/m) charge rms-radius is 
quite good. 

The magnetic form factor of ^Li has first been measured by Rand et al. |P45|| ; the 180° 
degree data, which reach 1.4/m~^, still had large statistical errors. Precise data were 



obtained by Lapikas et al. ||246|| at low g, and Bergstrom et al. [|247|| by an experiment 
performed using the 180° setup at Bates. The ^Li magnetic form factor calculated by 
Wiringa and Schiavilla, displayed in fig. |33|b, does not show as good an agreement with 
the data for momentum transfers q > lfm~^. This discrepancy with the data is not due 
to two-body currents, which have a relatively minor effect as ^Li is an isoscalar system; 
one therefore must assume that it is related to deficiencies in the VMC wave function. 
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Figure 34: Exp{S) calculation of the oxygen effective form factor (cr/crMott) 
by Mihaila and Heisenberg ||241|| , compared to the data of Sick and McCarthy 
p49| and Bentz [|5l. 



For Oxygen, Mihaila and Heisenberg ||241|| have recently used the exp{S) coupled- 
cluster expansion approach. The procedure is similar to the exp{S) method used at the 
time by Gari et al. [|14|] ( see fig. ^) who, in addition to ^ife, also considered heavier 
magic nuclei such as oxygen and calcium. In the coupled cluster approach, the nuclear 
wave function is expanded in the many-body Hilbert space in terms of multiconfigurational 
creation- and destruction operators. The expansion coefficients carry the information on 
the nuclear correlations, and no distinction between short- and long-range correlation 
needs to be made. The wave function is expanded in a 50huj space, and relative angular 
momenta up to 1=11 are taken into account. In this calculation, the Argonne V18 N-N 
interaction together with the Urbana IX three-body force is employed. Two-body currents 
are also included, and have a moderate effect at momentum transfers above 2 fm~^. The 
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Coulomb distortion is accounted for by Fourier-transforming the calculated form factor to 
r-space, and then calculating the Coulomb-distorted cross section (effective form factor) 
using a phase-shift code. 



The calculation is compared in fig. ^ to data |P49|| obtained in an experiment that 
used a windowless ice-target at LN2 temperature. The low-g data (0.3 - 0.5 fm~^) have 
been measured at Darmstadt by Bentz ||250|| . 

The agreement between calculation and experiment is amazingly good. In r-space, 
somewhat larger deviations are visible as the calculated rms-radius of 2.62/m [p51|| is 



smaller than the experimental one of 2.712±0.005/m. It will be interesting to see whether 
this exp{S) approach of Mihaila and Heisenberg can with similar success predict the 
properties of other p-shell nuclei. 



7 Conclusions and outlook 

In this paper we have reviewed the present status of the elastic form factors of nuclei with 
mass number A=2,3,4,6 and 16, both in terms of the cross section data available, the 
form factors extracted and their theoretical understanding. For the case of the deuteron, 
we have included the transition form factor to the singlet-S state. 

We have witnessed during the last years important progress in the theory used to 
calculate few-body wave functions and the corresponding electromagnetic form factors. 
The "nuclear standard model" — based on the assumption that nuclei are basically non- 
relativistic systems containing nucleons bound by the N-N interaction known from N-N 
scattering — is surprisingly successful. For the nuclei where the computational techniques 
are well under control, basically A < 4, the agreement between calculated ground state 
properties, and form factors in particular, is very satisfactory. 

The most important uncertainties for A < 4 still result from the poorly constrained 
nature of the N-N interaction off-shell. There, different N-N potentials give quite different 
answers. The ensuing differences mainly affect D-state properties. The same quantities 
are somewhat uncertain as the data base that fixes the S-D transition in the NN-system, 
represented by the ei-parameter, is still not very satisfactory and leaves too much free- 
dom for the N-N potentials. Here, only double-polarization experiments on n — p elastic 
scattering can be expected to lead to significant progress in the future. 

For the calculation of the electromagnetic form factors, it also has been very important 
to get the contributions of two-body terms, required by gauge invariance, under control. 
Much progress in this direction has been made, as many of the two-body terms now can 
be derived consistently with the N-N interaction employed. Diagrams not constrained by 
current conservation occasionally still involve significant uncertainties, though. 

One obvious avenue for a deeper understanding of light nuclei consist in a relativistic 
treatment. Given the typical momenta of nucleons in nuclei, a relativistic description is 
clearly desirable. A number of different approaches have been developed already, and 
we have addressed some of them when discussing the deuteron. For the heavier systems, 
the covariant calculations up to now have provided little in terms of form factors, the 
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quantities that yield the perhaps most detailed test. The covariant approaches also give 
us a better control over the off-shell behaviour and thus have the potential to partly 
do away with one ingredient of the non-relativistic calculations that still is somewhat 
unsatisfactory, the (phenomenological) three-body interaction. In this area, significant 
advances can be expected during the coming years. 

On the side of experiment, good progress has also been made. In particular, we now 
have a reasonably complete set of polarization data for electron-deuteron scattering that 
allows us to separate the deuteron monopole- and quadrupole form factors. For a long 
time, we had to contend with the unseparated structure function A{q), which is much less 
sensitive to the underlying physics. 

In this review, we have used for the nuclei A=2-4 the form factors resulting from a 
fit to the world data, such as to base the discussion on the most precise experimental 
information presently available. The numerical values of the various form factors, which 
partly have been determined in the course of writing this review, are available upon 
request. 

While our knowledge of the proton charge- and the neutron magnetic form factors 
has greatly improved during the last years, the neutron charge form factor Gen remains 
a continued source of uncertainty for the understanding of few-body charge form factors. 
The first double-polarization experiments now have produced results and lead to values 
of Gen that are more reliable than what was available up to now, but the uncertainties on 
Gen are still comparatively large, and the information is limited to moderate values of the 
momentum transfer. Here, important further progress is required in the future. With the 
CW-facilities today available and the progress in polarized target and polarimeter design, 
the necessary improvements are within reach. 

As far as the heavier nuclei go, it is very satisfactory to see the emergence of "exact" 
calculations (in the sense of the nuclear standard model) for p-shell nuclei. In the p-shell, 
these calculations can be used to "calibrate" the mean- field calculations that are standard 
(and often the only feasible ones) for heavier nuclei. 

When looking at the overall comparison of theory and experiment, a few notable 
discrepancies remain. Even the most modern calculations do not very well describe the 
^He magnetic form factor; the reason for this is not known and cannot easily be blamed on 
the two-body contributions which for this observable are largely of the model-independent 
kind. The comparison between the positions of the minima, and amplitudes of the form 
factor maxima, for the A— 2,3 systems also shows a persistent discrepancy. It is not clear 
if this can be taken as a signal for "new physics" beyond the nuclear standard model. It is 
worrisome as well that the calculations do not yet accurately explain the static properties 
{Q, n) of the simplest nucleus, the deuteron. 

It would evidently be of great interest to go to larger momentum transfer for some of 
the form factors. An extension of the ^He magnetic form factor to larger q will be carried 
out at JLAB, a significant extension for the two Tritium form factors is clearly feasible 
with available liquid-target technology. Pushing the T20 measurements for the deuteron 
to higher q will presumably require the development of a more efficient polarimeter of as 
yet unknown type. 
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Several aspects have received little attention in this review. Changes of the miclcon 
upon integration in the nucleus have not been discussed, essentially because there is little 
experimental evidence for it. We have also only occasionally discussed an understanding 
of the light nuclei in terms of quark constituents; the description in terms of nucleons is 
very successful, and clear signals for the need of a description in terms of the underlying 
quark degrees of freedom are hardly available. These signals probably only can be found 
by going to larger momentum transfers where the conventional description of light nuclei 
in terms of nucleonic constituents interacting via the N-N force must break down. 
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